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OVERVIEW 

Tensile  stress-strain  behavior  has  been  widely  recognized  as  the  most  important 
single  property  for  ceramic  matrix  composites.  In  this  study,  the  tensile  failure  mechanism 
of  the  ceramic  matrix  composites  has  been  investigated. 

First,  the  dog- bone  tensile  specimens  made  of  unidirectional  ceramic  matrix 
composite  namely:  Nicalon  Fiber  (S^C^,  (3  phase)  and  CAS  II  (Ca0-Al203-2Si02) 
matrix  with  fiber  volume  fractions  of  30%  and  40%  were  tested  in  room  as  well  as  in 
elevated  temperatures.  All  tests  were  conducted  inside  a  chamber  of  a  scanning  electron 
microscope  (SEM).  Therefore,  images  at  various  magnifications  revealing  the  intermediate 
tensile  damage  events  in  the  composite  specimen  were  obtained  in-situ  along  with  other 
testing  data.  Since  the  lens  in  the  SEM  is  always  kept  at  least  one  inch  away  from  the 
specimen  during  testing,  with  magnifications  ranging  from  50x  to  as  high  as  3000x,  the 
heat  on  the  specimen  won't  damage  the  electron  lens.  This  makes  it  possible  to  obtain 
images  while  doing  tensile  testing  in  high  temperature.  The  damage  histories  of  the 
specimen  at  different  temperatures  were  first  observed  on  the  monitor  of  the  SEM  and 
then  captured  immediately  either  by  Polaroid  camera  or  video  copy  processor.  Tensile 
stress-strain  curves  were  constructed  at  different  temperatures  (starting  from  room 
temperature,  then  250  °C,  400  °C,  600  °C  and  700  °C).  All  tests  exhibited  a  non-linear 
stress-strain  behavior.  The  effect  of  specimen  size  (meaning  specimens  with  different 
thickness  or  width  in  the  gauge  section)  on  the  failure  behavior  was  discussed.  It  was 
found  that  while  thin  specimens  (about  0.0625  "-0.07  "  in  either  thickness  or  width)  tend 
to  have  a  tail  at  failure,  retaining  some  load  carrying  capability,  thick  specimens  (greater 
than  0.0725"  on  both  thickness  and  width)  fail  catastrophically,  breaking  into  two  pieces. 
The  effect  of  temperature  on  tensile  failure  strength  was  found  not  quite  significant  within 
the  temperatures  tested.  Micrographs  taken  at  different  loading  and  temperatures  revealed 
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the  details  of  damage  process  from  first  matrix  cracking  to  multiple  matrix  cracking, 
sequential  random  fiber  slipping  and  fiber  breaking  to  eventual  composite  failure  by 
sudden  fiber  pull-out  and  collective  fiber  breaking  somewhere  in  the  gage  section.  Matrix 
crack  density  versus  the  stress  with  different  fiber  volume  fraction  and  at  different 
temperatures  was  also  studied.  The  results  indicate  that,  for  specimens  with  30%  fiber 
volume  fraction,  the  matrix  crack  initiation  stress  is  less  than  that  of  40%  fiber  volume 
fraction.  And  the  matrix  crack  density  (defined  as  the  number  of  matrix  cracks  per  mm) 
for  specimens  of  40%  fiber  volume  fraction  is  higher  than  that  for  those  with  30%  fiber 
volume  fraction.  It  is  also  noted  that  matrix  crack  initiation  stress  was  lower  in  high 
temperature  than  in  room  temperature  for  specimens  of  both  30%  and  40%  fiber  volume 
fraction. 

Then,  analytical  models  based  on  finite  element  method  and  the  singular  integral 
equation  technique  were  used  to  explain  the  tensile  damage  behavior  of  the  ceramic  matrix 
composites.  The  H-shaped  crack  geometry  was  used  in  the  finite  element  model  to  study 
the  interface  progression  and  its  effect  on  the  tensile  damage  behavior.  The  results  from 
finite  element  model  compare  well  with  the  experimental  results.  The  singular  behavior  at 
the  transverse  and  interface  crack  tips  were  studied  using  singular  integral  equations.  The 
stress  intensity  factors  and  strain  energy  release  rates  were  calculated  for  various  crack 
geometry  and  were  used  to  explain  the  failure  mechanism  of  the  composites.  The  results 
from  the  singular  integral  equation  technique  predict  that  once  the  transverse  matrix 
cracks  are  formed,  they  will  propagate  to  the  fiber/matrix  interface.  This  behavior 
conforms  to  the  observed  behavior.  Both  models  assume  that  the  composite  consists  of 
equal  spaced  fiber  strips  in  the  matrix  material  and  the  problem  is  simplified  as  two 
dimensional. 


3 


1  INTRODUCTION 

The  research  on  ceramic  matrix  composite  materials  has  intensified  in  recent  years 
due  to  some  of  the  appealing  features  of  ceramics;  such  as  great  stability  and  resistance  to 
oxidation  under  hostile  (high  temperature  or  corrosive)  environments.  In  comparison  to 
their  metallic  and  polymer  counterparts,  ceramics,  being  brittle  and  low  in  tensile  strength 
and  fracture  toughness,  traditionally  have  had  little  use  in  structural  applications. 
However,  when  reinforced  with  fibers,  ceramic  matrix  composites  exhibit  an  increase  in 
fracture  toughness  and  tensile  strength  in  room  as  well  as  high  temperatures  [1-2].  A 
variety  of  ceramic  matrix  composite  systems  have  been  or  are  being  developed  for 
engineering  applications  ranging  from  cutting  tools  to  aerospace  structures  [3-4].  For 
examples,  one  can  find  studies  on  the  following  systems  of  ceramic  matrix  composites: 
C/glass  [5],  C/SiC  [6],  SiC/glass  [7],  SiC/LAS  glass  ceramic  [8-13],  SiC/BMAS  glass- 
ceramic  [13],  SiC/alumina  [14],  SiC/mullite  [15],  SiC/SiC  [16]. 

In  the  aforementioned  studies,  experimental  results  on  damage  behaviors  of  ceramic 
matrix  composites  were  obtained  by  failing  the  specimens  with  either  tensile  or  three-point 
bending  loading  at  room  or  elevated  temperature  [9-12].  During  the  thermomechanical 
testing,  the  load  vs.  the  displacement  (and  hence  the  stress  vs.  the  strain)  curves  were 
recorded  but  only  the  postmortem  damage  patterns  were  identified  by  either  a  scanning 
electron  or  an  optical  microscope.  It  is  well  known  that  the  stress-strain  relations  of  most 
ceramic  matrix  composites  under  thermomechanical  loading  usually  exhibit  nonlinear 
behaviors.  For  such  ceramic-matrix  composites,  the  failure  mechanisms  that  cause  this 
nonlinear  stress-strain  relationship  is  more  complex  than  that  of  their  monolithic 
counterparts.  They  come  from  the  results  of  multiple  matrix  cracking  and  sequential  fiber 
breaking  due  to  weak  interfacial  bonding  between  the  fibers  and  the  matrix.  The 
approaches  described  above  although  proved  important  data  to  the  overall  understanding 
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of  the  failure  mechanisms  of  the  ceramic  matrix  composites,  they  fail  to  correlate  the 
nonlinear  stress-strain  (and  hence  the  stiffness  reduction)  behavior  with  the  intermediate 
damage  progression  events  inside  the  ceramic  matrix  composite  specimen.  A  good 
experimental  approach  should  provide  information  on  the  damage  history  of  the  specimen 
that  can  be  used  to  correlate  the  nonlinear  stress-strain  (and  hence  the  stiffness  reduction) 
behavior  with  the  intermediate  causative  damage  events  that  occurred  on  the  ceramic 
matrix  composite  specimen.  Recently,  in  [17-18],  the  tensile  behavior  of  a  Nicalon/CAS  II 
system  was  studied.  Damage  patterns  were  identified  and  micrographs  were  taken  to 
capture  the  matrix  crack  propagation.  However,  these  micrographs  were  not  taken  at  the 
same  location  and  therefore  cannot  truly  correlate  the  damage  progression  with  the 
nonlinear  stress-strain  behavior.  And  those  studies  dealt  only  with  room  temperature.  In 
this  study,  experiments  were  conducted  inside  the  chamber  of  a  scanning  electron 
microscope  equipped  with  a  custom  designed  tensile/heating  substage.  This  made  it 
possible  to  directly  observe  and  record  in  situ  the  progressive  tensile  damage  behavior  of 
the  ceramic  matrix  composites  from  the  very  first  matrix  crack  to  complete  fracture  of  the 
specimen  at  any  location  in  the  gage  section.  Test  results  for  both  room  and  elevated 
temperatures  were  obtained.  One  particular  advantage  of  using  SEM  is  for  high 
temperature  testing.  Since  in  SEM,  unlike  in  optical  microscope,  the  electronic  lens  is 
always  kept  at  least  one  inch  away  from  the  heated  specimen,  images  can  be  obtained  at 
magnifications  as  high  as  3000x.  without  damaging  the  electronic  lens  during  high 
temperature  testing.  This  makes  the  technique  very  appealing  for  high  temperature  testing. 

Among  the  factors  that  affect  the  overall  strength  and  toughness  of  the  ceramic 
matrix  composites,  the  following  are  generally  considered  the  most  important:  (a)  the 
thermoelastic  properties  of  the  constituents  (i.e.,  the  matrix  and  the  fibers);  (b)  the  relative 
strength  of  the  interface  between  fibers  and  matrix;  (c)  the  volume  fraction  and 
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airangement  of  the  fibers;  and  (d)  the  ambient  temperature.  Based  on  the  above 
considerations  and  observations  during  the  test,  two  analytical  models  have  been 
developed.  One  model  utilizes  the  finite  element  technique  and  is  based  on  a  periodic 
multiple-row  H-shaped  crack  configuration.  The  other  uses  the  singular  integral  equation 
method  with  periodic  transverse  and  interface  crack  configuration.  In  both  models  the  real 
problem  is  formulated  in  two  dimensional  domain. 

2  THE  EXPERIMENTAL  WORK 

The  experimental  work  is  vital  in  this  study.  It  provides  important  insight  and  data 
for  the  understanding  of  the  failure  process  of  this  type  of  composite  material.  And  it  also 
provides  the  foundation  upon  which  analytical  models  are  based.  The  technique  used  in 
this  study  is  believed  to  be  first  developed  in  the  area  of  experimental  mechanics  research. 

2.1  THE  TESTING  PROCEDURE 

Figure  2.1  shows  a  schematic  drawing  of  the  experimental  set-up  which  uses  a 
scanning  electron  microscope  equipped  with  a  tensile  and  heating  substage  to  perform  the 
micromechanical  tensile  testing  of  a  ceramic  matrix  composite  specimen  under  high 
temperature.  In  this  study,  a  Hitachi  S-2400  scanning  electron  microscope  is  used  which 
is  equipped  with  a  custom  designed  E.F.  Fullam  tensile/heating  combined  substage.  The 
specimen  used  is  simple  “dog-bone”  shaped  specimen  such  as  the  one  shown  in  Fig.  2.2. 
The  actual  experimental  set-up  and  tensile/heating  substage  are  shown  in  Fig.  2.3.  The 
material  used  in  this  study  (Nicalon-fiber/CAS  D  matrix  composite)  was  obtained  from 
Coming  Glass  Works.  Table  2.1  shows  the  thermomechanical  properties  of  the 
constituents  of  the  composite.  Specimens  of  both  30%  and  40%  fiber  volume  fractions 
were  tested.  To  study  the  size  effect,  we  varied  either  the  thickness  or  the  width  of  the 
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specimen  for  specimens  with  40%  fiber  volume  fraction.  Fibers  were  unidirectionally 
aligned  with  the  gage-length  direction.  The  tensile  specimens  were  made  by  first  cutting 
the  large  rectangular  ceramic  matrix  composite  panel  as  supplied  by  Coming  Glass  Works 
into  smaller  rectangular  plates  using  a  Leco  VC-50  Cari/Cut  fine-mesh  diamond  saw. 
Then  the  gage  section  of  the  specimens  was  shaped  using  the  same  diamond  saw  with 
special  a  holder.  The  gage  section  of  the  specimens  was  grounded  to  its  dimensions  using 
a  Dremel  motorized  hand-held  grinder  with  silicon  carbide  and  alumina  oxide  stones. 


Table  1.  Properties  of  Nicalon  fiber  and  CAS  II  matrix 


To  conduct  high  temperature  testing,  the  bottom  surface  of  the  central  part  of  the 
specimen  as  shown  in  Fig.  2.1,  was  placed  in  direct  contact  with  the  E.F.  Fullam  heating 
element  (called  heater)  which  is  a  rectangular  plate  made  of  ceramic  material  with  fuse 
wire  circuit  inside.  The  maximum  heated  area  of  the  heating  element  is  0.65"x0.25".  It 
can  sustain  a  maximum  working  temperature  of  1100°C  and  is  equipped  with  a  water- 
cooled  heat  sink  for  continuous  operation.  Temperature  is  measured  by  three  platinum 
30%  rhodium-platinum  6%  rhodium  thermocouples  and  controlled  by  a  stand-along  DC 
power  supply  with  adjustable  voltage  and  current  knobs.  Since  the  whole  operation  is 
conducted  inside  the  chamber  of  a  scanning  electron  microscope  (in  this  study  a  Hitachi  S- 
2400)  which  is  usually  vacuumed  at  1.5xl0‘6  Pa  or  better,  no  heat-loss  will  occur  due  to 
thermal  convection.  The  top  surface  of  the  specimen  was  first  polished  using  fine-grade 
diamond  paste  containing  15-,  6-,  and  1-micron  particles  until  the  surface  was  well 
finished  and  the  fibers  and  matrix  could  be  seen  clearly  under  a  Nikon  UM-2  microscope. 
Then  the  specimen  was  cleaned  using  a  Bronson  ultrasonic  cleaner  for  10  minutes.  Finally 
the  bottom  surface  of  the  specimen  was  coated  with  silver  paint  to  prevent  electric 
charging  and  to  achieve  better  image  before  it  was  mounted  into  the  E.F.  Fullam  tensile 
substage.  Because  it  is  very  difficult  to  drill  holes  in  a  ceramic  composite  specimen,  the 
top  and  bottom  ends  of  the  specimen  were  mounted  through  stainless  steel  clamps  with 
serrated  teeth  to  the  crossheads  of  the  E.F.  Fullam  tensile  substage.  To  help  in  alignment 
and  prevent  slippage  during  testing  between  the  ceramic  composite  specimen  and  the 
serrated  clamps,  cyanoacrylate-based  extra- strength  epoxy  was  also  applied  on  the  clamp- 
specimen  interfaces.  The  clamped  specimen  was  then  mounted  to  the  tensile  substage  and 
was  left  to  cure  for  at  least  20  hours  to  ensure  that  the  epoxy  had  hardened  and  was 
completely  dry  before  testing. 
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Finally,  the  specimen  together  with  the  tensile/heating  assembly  was  placed  inside 
the  chamber  of  a  Hitachi  S-2400  scanning  electron  microscope  which  is  also  equipped 
with  a  backscatter  detector  to  enhance  the  image.  By  using  the  X-Y  staging  control  of  the 
scanning  electron  microscope,  micrographical  patterns  within  the  central  gage  area  can  be 
observed  and  recorded.  The  tensile  substage  is  driven  by  a  variable- speed  motor  which 
has  a  maximum  speed  of  90  rpm.  Through  a  gear  mechanism  of  100: 1  reduction  ratio,  the 
crosshead  speed  can  be  controlled  within  0.127  mm/min  (or  0.005  in/min).  The  tensile 
stage  can  provide  a  tensile  load  of  up  to  455  kgs  (1000  lbs).  The  applied  load  was 
increased  gradually  until  the  specimen  failed  totally.  The  damage  to  the  composite  as  the 
load  increased  was  first  observed  on  the  monitor  of  the  SEM.  Then  a  sequence  of 
micrographs  were  taken  to  get  hard  copies  of  various  microcracking  and  damage  patterns 
of  the  specimen  at  different  loading  levels.  Since  the  specimen  is  very  thin  thickness 
(ranges  from  1.0  mm  to  2.0  mm  or  0.04"  to  .08"),  it  is  expected  that  failure  will  occur 
through  the  thickness.  Thus  the  recorded  micrographs  of  the  surface  fracture  can 
represent  through-the-thickness  failure  of  the  specimen.  Quality  of  the  image  and  hence 
the  quality  of  the  micrographs  was  further  enhanced  by  transmitting  the  image  signal  to  a 
computer  which  is  equipped  with  an  Imaging  Technology  Advanced  Frame  Grabber 
(AFG)  digital  image  analyzer  hardware.  With  the  help  of  the  installed  software,  sharper 
images  were  obtained  through  the  contrast  and  edge  enhancement  operations.  And  text 
can  be  added  to  the  micrographs. 

The  tensile  stage  is  designed  in  such  a  way  that  when  the  load  increases,  the  top  and 
bottom  crossheads  move  in  opposite  directions  to  minimize  the  shift  of  the  observed  site. 
This  is  achieved  by  machining  the  stainless-steel  loading  columns  into  worms  of  reverse 
directions.  Thus  searching  and  refocusing  the  damaged  zone  after  each  load  increment  are 
very  handy.  The  applied  loads  were  recorded  by  a  miniature  load  cell  equipped  with  a 
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digital  readout.  The  relative  displacements  of  the  crossheads  were  measured  by  a  high- 
precision,  strain-gage-type  extensometer  which  is  also  fitted  with  a  digital  readout.  The 
load-displacement  data  were  recorded  and  converted  into  a  stress-strain  curve.  The 
evolution  of  damage  recorded  by  micrographs  were  identified  with  the  corresponding 
stress  and  strain.  Finally  the  ruptured  specimens  were  observed  under  a  Nikon  UM-2 
universal  measuring  microscope  for  further  postmortem  examination. 

2.2  EXPERIMENTAL  RESULTS 
2.2.1  Stress-Strain  Relations 

Typical  stress-strain  curves  for  the  room-temperature,  tensile  damage  behavior  of  a 
Nicalon/CAS  II  with  30%  and  40%  fiber  volume  fractions  are  shown  in  Fig.  2.4  and  Fig 
2.5  respectively.  As  depicted  in  these  two  figures,  the  tensile  damage  behavior  of  the 
ceramic  matrix  composite  specimen  is  characterized  by  a  nonlinear  curve  made  up  of 
three  sections.  The  characteristics  are  the  same  for  all  tests  conducted  in  this  study,  even 
though  microscopically  these  specimens  might  look  quite  different.  Upon  load  application, 
the  relation  between  stress  and  strain  was  linear  and  its  slope  was  equivalent  to  the 
stiffness  of  an  intact  Nicalon/CAS  II  specimen  (18.15  Msi  or  125  GPa  for  Vf  =  30%  and 
19.6  Msi  or  135  GPa  for  Vf  =  40%  ).  For  specimens  with  30%  fiber  volume  fraction,  the 
slope  changed  at  about  25  ksi  (Point  A  in  Fig  2.4).  While  for  40%  fiber  volume  fraction, 
the  slope  changed  at  about  30  ksi  (Point  A  in  Fig.  2.5).  First  matrix  crack  is  believed  to 
start  at  or  slightly  below  of  point  A  in  both  Vf  =  30%  and  40%  cases.  This  can  be  further 
inferred  from  the  relationship  between  matrix  crack  density  and  tensile  stress  as  will  be 
discussed  later  in  this  section.  In  most  cases,  matrix  cracking  initiates  either  at  the  edge  or 
at  a  location  where  the  fiber  spacing  is  maximum.  Fig  2.6  shows  some  micrographs  taken 
during  one  test.  The  second  micrograph  in  Fig.  2.6  clearly  indicates  that  matrix  crack 
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initiated  not  from  the  voids  but  from  the  location  where  spacing  between  the  two  adjacent 
fibers  is  the  largest.  While  in  Fig.  2.7,  matrix  cracks  initiated  from  both  the  edge  and  a 
location  where  fiber  spacing  is  very  large  if  not  the  largest.  The  matrix  cracks  then 
propagate  perpendicular  to  the  fiber  direction  throughout  the  whole  width  of  the  gage 
section  to  form  multiple  matrix  cracks.  Unlike  the  monolithic  materials  where  failure  is 
controlled  by  a  critical  crack  size,  the  failure  of  the  composite  goes  through  a  process  of 
damage  accumulation.  The  composite  is  insensitive  to  the  voids  and  can  tolerate  very  large 
cracks  before  failure.  At  point  B  (35  ksi  for  Vf  =  30%  and  40  ksi  for  Vf  =  40%),  the 
development  of  multiple  matrix  cracks  reached  a  saturated  stage.  Regularly  spaced  matrix 
cracks  were  formed  in  the  whole  gage  section  of  the  specimen.  This  happened  with  only  a 
small  increment  of  tensile  stress  (about  10  ksi,  as  can  be  seen  from  point  A  to  point  B  in 
both  Fig.  2.4  and  Fig.  2.5).  Further  increase  of  tensile  stress  creates  no  or  very  little 
additional  matrix  cracking.  Matrix  crack  opening,  fiber  debonding,  breaking  and  slipping 
will  dominate  the  rest  of  the  failure  process.  Upon  reaching  point  C  ( about  62  ksi  for  Vf  = 
30%  and  58  ksi  for  Vf  =  40%  ),  one  surface  of  the  matrix  cracks  in  the  gage  area  started 
to  open  up  with  crack  opening  displacement  far  more  larger  than  the  rest  of  the  matrix 
cracks.  This  then  was  accompanied  by  massive  fiber  breaking  and  pull-out  in  that  surface. 
And  the  eventual  separation  of  the  specimen  caused  the  load  to  drop  substantially. 
Depending  on  whether  the  specimen  is  thin  (either  in  thickness  or  in  width)  or  thick,  there 
might  be  a  tail  in  the  stress-strain  curve  at  the  load  drop.  Thin  specimens  tend  to  have  a 
tail  at  failure  as  indicated  in  Fig.  2.5  where  both  thin  and  thick  specimens  were  tested.  One 
explanation  is  that  thin  specimens  are  prone  to  bending  during  the  test.  This  might  be 
responsible  for  the  lower  failure  strength  and  a  tail  at  failure  for  thin  specimens.  For  thick 
specimens,  failure  is  always  catastrophic.  Fig.  2.8  shows  the  progression  of  damage  for  a 
thin  specimen.  At  failure,  the  specimen  (this  occurred  only  for  thin  specimens)  was  kept  in 
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one  piece  by  the  fibers.  Further  loading  after  failure  was  possible  as  shown  in  Fig.  2.9. 
Micrographs  in  Fig.  2.10  however,  show  the  catastrophic  failure  of  a  thick  specimen. 
Upon  reaching  the  maximum  stress  value,  the  specimen  fails  catastrophically  with  a  very 
big  crack  opening  at  the  failure  surface.  Observation  on  the  failure  surfaces  of  the 
specimens  with  30%  and  40%  fiber  volume  fractions,  indicated  that  there  are  significant 
differences  in  both  the  amount  and  the  length  of  fiber  pull-out  between  the  two  specimens. 
For  specimens  with  Vf  =  30%,  both  the  amount  and  the  length  of  fiber  pull-out  at  failure 
surface  are  less  than  those  with  Vf  =  40%  as  can  be  seen  from  Fig.  2.1 1  and  Fig.  2.12.  The 
lesser  amount  of  fiber  pull-out  and  the  shorter  fiber  pull-out  length  is  an  indication  of 
stronger  interfacial  bonding  strength.  This  explained  why  the  tensile  failure  strength  of  Vf 
=  30%  is  slightly  higher  than  that  of  Vf  =  40%.  The  smoothness  of  fiber  pull-out  surface 
as  shown  in  Fig.  2.13  is  evidence  of  non-chemical  bonding  between  fibers  and  matrix 
which  generally  implies  weak  interfacial  bonding  strength.  The  fiber  pull-out  length  and  its 
amount  also  deviated  substantially  among  the  same  batch  of  specimens  as  shown  in  Fig. 
2.14,  implying  that  the  fiber/matrix  interfacial  strength  may  vary  for  the  same  batch  of 
ceramic  matrix  composite  specimens.  This  is  also  reflected  by  the  fluctuation  of  the 
maximum  tensile  stresses  (Point  C)  among  these  tests. 

Figure  2.15  is  a  set  of  micrographs  showing  typical  damage  patterns  of  the  specimen 
at  different  magnifications  after  the  stress  reached  point  B  in  the  stress-strain  curve.  At 
close-up,  one  can  clearly  identify  the  matrix  crack  opening,  fiber  breaking,  slipping  and 
fiber  bridging  of  the  matrix  cracks.  Also  from  Fig  2.15,  one  can  observe  the  typical  "H" 
shaped  crack  pattern  formed  by  intersection  of  the  transverse  matrix  cracks  with  interface 
cracks  during  the  failure  process.  This  "H”  shaped  crack  configuration  will  be  used  in  the 
analytical  models  discussed  later.  Besides  the  easily  identified  matrix  crack  patterns,  fiber 
breaking  also  exhibits  some  patterns  as  shown  in  Fig.  2.16  and  Fig.  2.17.  In  Fig.  2.16,  one 
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fiber  breaks  at  a  location  away  from  the  matrix  cracks  first,  then  a  ray  of  breaking  fibers  is 
form  along  a  line  slanted  away  from  the  line  of  matrix  cracks.  Fiber  breaking  occurred 
mosdy  after  the  stress  reached  point  B,  i.e.  after  the  multiple  matrix  cracks  have  been 
formed.  Micrograph  in  Fig  2.17  gives  another  pattern  of  fiber  breaking  in  which  fiber 
breaking  is  randomly  distributed  throughout  the  whole  gage  section. 

The  stress-strain  relations  of  the  Nicalon/CAS  II  composite  at  high  temperatures 
will  be  discussed  later. 

2.2.2  Matrix  Crack  Density 

Matrix  crack  density  defined  as  the  number  of  matrix  cracks  per  1  mm  length  in 
fiber  direction,  was  used  to  characterize  the  failure  process  of  the  composite.  With  the 
help  of  the  scanning  electron  microscope  equipped  with  a  backscatter  detector,  images 
reflecting  the  tensile  damage  pattern  of  the  composite  from  the  first  matrix  crack  to  the 
eventual  failure  of  the  composite  were  captured.  Fig  2.18  and  Fig.  2.19  are  some 
micrographs  showing  the  development  of  matrix  cracks  with  increasing  stresses  for  Vf  = 
30%  and  Vf  =  40%  respectively.  Based  on  the  number  of  cracks  counted  in  the  frame 
shown  in  these  micrographs  and  the  associated  stress  level,  the  matrix  crack  density  versus 
tensile  stress  curve  was  constructed  as  shown  in  Fig  2.20  and  Fig  2.21  for  30%  and  40% 
fiber  volume  fractions  respectively.  From  the  matrix  crack  density  curves,  it  is  seen  that 
the  matrix  crack  initiation  stress  is  25  ksi  for  Vf  =  30%  and  30  ksi  for  Vf  =  40%.  These 
results  happen  to  coincide  with  point  A  on  the  stress- strain  curves  as  shown  in  both  Fig 
2.4  and  Fig  2.5. 


2.2.3  Effect  of  Fiber  Volume  Fractions 
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As  can  be  seen  from  the  stress-strain  and  matrix  crack  density  curves  for 
specimens  of  both  30%  and  40%  fiber  volume  fractions,  the  fiber  volume  fraction  affects 
the  tensile  behavior  substantially.  Matrix  crack  initiation  stress  is  less  for  specimens  with 
less  fiber  volume  fraction  (compare  points  A  in  Fig  2.4  and  Fig.  2.5).  This  result  is  in 
agreement  with  that  given  in  [19].  But  the  tensile  failure  strength  for  Vf  =  30%  is  higher 
than  that  for  Vf  =  40%  (points  C  in  the  stress-strain  curves).  This  might  be  the  result  of  a 
relatively  stronger  interfacial  bonding  strength  for  specimens  with  Vf  =  30%.  The  matrix 
crack  density  was  found  to  increase  with  increasing  fiber  volume  fraction  of  the  specimen. 
This  finding  also  agrees  with  that  reported  in  [19]. 

2.2.4  Temperature  Effects 

The  same  specimens  were  tested  in  the  SEM  at  250  °C,  400  °C,  600  °C  and  700 
°C.  Fig.  2.22  and  Fig.  2.23  show  the  stress-strain  curves  obtained  at  higher  temperatures 
for  fiber  volume  fractions  30%  and  40%  respectively.  Within  the  temperature  range  tested 
(  room  temperature,  250  °C,  400  °C,  600  °C  and  700  °C),  no  significant  changes  in 
ultimate  tensile  strength  have  been  observed  in  the  stress-strain  curves.  Fig.  2.24  and  Fig. 
2.25  show  the  temperature  effects  on  the  ultimate  tensile  failure  strength  of  the  specimens 
for  both  30%  and  40%  fiber  volume  fractions.  There  was  a  slight  increase  in  tensile  failure 
strength  for  Vf  =  30%  at  T  =  250  °C  .  Micrographs  showing  the  final  failure  surfaces  of 
the  specimens  with  Vf  =  30%  at  different  temperatures  are  compared  in  Fig.  2.26.  It  is 
seen  that  both  the  amount  and  the  length  of  fiber  full-out  are  the  smallest  for  T  =  250  °C. 
This  might  explain  why  the  tensile  strength  is  maximum  at  T  =  250  °C.  For  Vf  =  40%,  this 
increase  of  tensile  strength  happened  at  T  =  400  °C.  The  phenomenon  of  slight  increase  of 
tensile  strength  might  be  the  result  of  release  of  thermal  residual  stress  and  changing  of 
interfacial  bonding  strength  due  to  temperature.  Matrix  crack  density  did  not  show 
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significant  changes  with  temperatures  for  Vf  =  30%.  Micrographs  in  Fig.  2.27  show  the 
development  of  matrix  cracks  at  T  =  600  °C.  Fig.  2.28  shows  the  relation  between  the 
matrix  crack  density  and  temperatures  (ranging  from  room  to  600  °C)  for  Vf  =  30%.  But 
for  Vf  =  40%,  it  was  found  that  matrix  crack  density  decreased  with  increasing 
temperature  as  shown  in  micrographs  of  Fig.  2.29  taken  at  T  =  400  °C.  Fig.  2.30  plots  the 
relation  of  matrix  crack  density  versus  temperature  for  Vf  =  40%.  Matrix  crack  initiation 
stress  also  decreased  slighdy  at  high  temperatures  for  Vf  =  40%.  Micrographs  in  Fig  2.31 
and  Fig.  2.32  show  the  fiber  breaking  and  slipping  process  in  high  temperatures.  Since  this 
composite  is  basically  a  glass  based  material,  the  stress-strain  behavior  is  expected  to  be 
quite  different  when  temperature  reaches  over  800  °C. 

2.2.5  Specimen  Sizing  Effects 

To  study  the  size  effect,  specimens  of  40%  fiber  volume  fraction  with  different 
thickness  and  width  combination  were  used  in  the  tests.  It  was  found  that  specimens  with 
either  relatively  thin  thickness  or  width  (about  0.0625-0.07")  tended  to  have  a  quite 
different  failure  pattern  than  that  with  thicker  one  (greater  than  0.0725"  in  both  thickness 
and  width).  For  thin  specimens,  first,  the  tensile  strength  was  quite  lower  than  that  of  thick 
specimens;  second,  when  the  specimen  failed,  it  had  a  tail  in  the  stress-strain  curve 
indicating  that  it  still  retained  some  load  carrying  capability  and  the  specimen  was  kept  in 
one  piece  by  the  fibers.  The  thick  specimens  on  the  other  hand  always  failed 
catastrophically.  A  comparison  of  tensile  strength  between  thin  and  thick  specimens 
against  temperature  is  shown  in  Fig  2.33.  The  discrepancy  between  the  two  cases  might  be 
the  result  of  bending  that  might  have  occurred  during  the  test  as  explained  earlier.  Size 
effect  of  the  specimens  on  the  flexural  strength  of  unidirectional  carbon  epoxy  composites 
has  been  reported  in  the  literature  [20-21].  Since  there  are  some  definite  relations  between 
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flexural  strength  and  tensile  strength,  tensile  strength  will  be  affected  by  the  size  of  the 
specimen  as  well.  From  the  standpoint  of  Weibull  failure  theory,  thick  specimens  tend  to 
contain  large  defects,  therefore  they  should  fail  at  a  lower  tensile  loading.  Since  the 
composite  with  a  relatively  weak  interfacial  bonding  strength  is  insensitive  to  the  defects, 
the  opposite  results  were  obtained  here.  The  specimen  didn't  fail  by  one  defect  or  crack 
when  it  reaches  the  critical  size,  but  through  the  damage  accumulation  all  over  the  gage 
section.  For  thin  Nicalon/CAS  II  specimens  with  40%  fiver  volume  fraction,  the  stress- 
strain  curves  at  high  temperatures  were  shown  in  Fig.  2.34. 
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3  THEORETICAL  MODELS 

As  explained  earlier,  during  testing,  it  was  found  that,  after  the  initiation  of 
transverse  matrix  cracks  at  point  A  (in  both  Figs.  2.4  and  2.5),  with  increasing  load, 
regularly  spaced  multiple  transverse  cracks  were  formed.  When  the  tensile  stress  reached 
a  certain  value,  almost  no  new  transverse  cracks  were  generated  until  failure.  During  this 
process  further  transverse  matrix  cracks  opening,  fiber  debonding,  breaking  and  sliding 
at  the  fiber/matrix  interface  is  believed  to  dominate  the  failure  process.  This  typical 
failure  feature  is  best  described  by  the  periodic  H-cracks  configuration  shown  in  Fig. 
2.15.  Two  models,  one  using  singular  integral  equation  technique,  and  the  other  using 
finite  element  method,  are  adopted  to  either  explain  the  failure  mechanism  or  simulate 
the  observed  tensile  behavior  of  ceramic  matrix  composites. 

The  problem  of  concern  is  basically  a  three  dimensional  problem  because  of  the 
distribution  of  fibers  in  the  matrix.  However,  if  one  assumes  that  the  fibers  are  made  of 
composite  strips  of  width  2HX,  the  problem  can  be  treated  in  two  dimensions. 

3.1  Singular  Integral  Equations  Formulation. 

It  has  been  shown  that  the  singular  integral  equation  technique  is  a  very  powerful 
tool  in  dealing  with  crack  problems.  Fig.  3.1  is  a  sketch  of  the  proposed  mathematical 
model.  It  is  assumed  that  fibers  are  equally  spaced  in  a  ceramic  matrix,  and  that  the 
thermomechanical  and  fracture  properties  of  the  fiber  and  the  matrix  are  known.  The 
model  contains  cracks  perpendicular  to  as  well  as  parallel  to  the  fibers  at  the  fiber/matrix 
interface  and  which  are  assumed  to  be  periodic.  By  choosing  different  geometrical 
parameters,  one  can  generate  various  crack  geometries.  For  example:  If  one  sets  ax=0, 
a2=H2  and  bx=0  in  Fig.  3.1,  then  the  model  reduces  to  a  geometry  with  periodic  H-shaped 
cracks  as  shown  in  Fig.  3.2.  This  configuration  closely  resembles  the  actual  cracked 
geometry  observed  in  testing. 
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3.1.1  Formulation  of  the  problem 

3.1.1.1  Equilibrium  equations 

The  equilibrium  equations  expressed  in  terms  of  displacements  for  orthotropic  materials 
are  as  follow: 


32u  32u  32v 

^3x2  +  3y2  +  ^33x3y 

(3.1a) 

32v  32v  32u 

3x2  +  ^3y2  +  ^33x3y 

(3.1b) 

where  pl5  p2  and  (33  are  material  constants  [23] 

Assume  that  the  solutions  are  of  the  following  forms: 

u(x,y)  =  ud)(x,y)  +  u<2>(x,y) 

(3.2a) 

v(x,y)  =  vd)(x,y)  +  v(2>(x,y) 

(3.2b) 

where 

2  r°° 

u(*)(x,y)  =  —  f(a,x)cosayda 

It  Jo 

(3.3a) 

2  r“ 

v(1)(x,y)  =  ~  g(a,x)sinayda 

It  Jo 

(3.3b) 

and 

2  r°° 

u(2)(x,y)  =  —  h(a,y)sinaxda 

It  Jo 

(3.4a) 

2  f°° 

v(2)(x,y)  =  —  {( a,y)cosaxda 

It  Jo 

(3.4b) 

here  f(a,x),  g(a,x),  h(a,y)  and  l( a,y)  are  shape  functions 

Substituting  (3.3ab)  and  (3.4ab)  into  (3. lab),  one  gets  the  following  characteristic 

equation: 

r4  +  P4r2  +  p52  =  0 

(3.5) 

where 


and 


P5  =VlVPi 
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and  the  roots  are  obtained  as: 

ri  =  -  r3  =  wj  +  iw2  =  "V^P 4+  Pg)/2 
r2  =  -  r4  =  w3  +iw4  =  V(-p4  -  pg)/2 
where  Pg  =  ^/p42  -4p52 


Then  we  have: 


f(a,x)  =  A(a)eriax  +  B(a)e'riax  +C(a)er2ax  +  D(a)e-r2«x 

g(a,x)  =  p7[A(a)eriax  -  B(oc)e‘riax] 

(3.6a) 

+  P8[C(a)er2ax  -  D(a)e'r2ax] 
h(a,y)  =  E(a)erlay/P5  +  F(a)e‘rlay/P5 

(3.6b) 

+  G(a)er2ay/P5  +  H(a)e-r2«y/p5 

Ka,y)  =  p9[E(a)eriay/P5  -  F(a)e'rlay/P5] 

(3.6c) 

+  p10[G(a)  er2«y/p5  -  H(a)e‘r2ay/p5] 

(3.6d) 

Here  two  types  of  material  will  be  distinguished  according  to  whether  the  roots  of 

the  characteristic  equation  are  real  or  complex. 

Material  type  I:  where  both  q  and  r2  are  real  numbers(w2=w4=0). 

Material  type  II:  where  q  and  r2  are  complex  numbers. 

Applying  symmetry  conditions: 

u(x,y)  =  -u(-x,y) 

(3.7a) 

v(x,y)  =  -v(x,-y) 

we  obtain: 

(3.7b) 

B(a)  =  -A(a)  D(a)  =  -C(a) 

(3.8ab) 

F(a)  =  -E(a)  H(a)  =  -G(a) 

(3.8cd) 

Material  type  I  will  be  considered  because  most  materials  fall  into  this  category.. 
Based  on  the  fact  that  both  u  and  v  vanish  when  y  00 ,  and  that  the  problem  of  concern 
is  symmetrical  about  the  x  axis,  the  displacement  functions  can  be  written  as: 
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4  r 

u(x,y)  =  —  I  [A(a)sinh(wi<xx)  +  C(a)sinh(w3ax)]cosayda 
71  Jo 

+  -  f"  [E(a)e-lwllay/P5  +  G(a)e-Iw3lay/P5]sinaxda 

K  Jo 

4 

v(x,y)  =  —  [P7A(a)cosh(wiax)  +  P8C(a)cosh(w3ax)]sinayda 
K  Jo 

--  p  [(39sign(w1)E(a)e-|wllay/p5  +  p10sign(w3)G(a)e-|w3lay/p5]cosaxda 
7C  Jo 


(3.9a) 


(3.9b) 


3.1.1.2  Stress-strain  relations 


Hxiiy 


then  for  orthotropic  materials,  under  plane  stress  conditions' 


°>=(x'y)=i^AEx+i^Aey 

(3.10a) 

^  v  V^v  1 

Oy(x,y)  —  g  ^tx  +  g  Aey 

X  X 

(3.10b) 

v; 

X 

•< 
s _ . 

II 

> 

(3.10c) 

where  Ex,  Ey,  Gxy,  vxy  and  Vyx  are  the  othortropic  material  constants. 

3.1.1.3  Strain-displacement  relations 


In  terms  of  u  and  v,  the  strains  are  given  as: 

3u  dv  ,  3u  5v 

£x  =  ax’  ey  “  ay  and  Yxy ~  (ay  +  ax} 


(3.11abc) 

Substituting  eqns  (3.9ab)  into  (3.11)  and  then  into  (3.10),  the  stress  can  be  written 


as: 


71(1  ^yVyx-gx(x,y)  =  [“^(aJe-'wi'ay/Ps  +  y2G(a)e-|w3lay/p5]acosaxda 
^x  Jo 

+  f  [2y3A(a)cosh(w1ax)  +  2y4C(a)cosh(w3ax)]acosayda  (3.12a) 
Jo 

7t(1-vXyvyX)gy(x,y)  =  J“[y5E(a)e-lw1lay/p5  +  y6G(a)e-Iw3lay/p5]acosaxda 
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+  f  [lyyACcOcoshCwiOtx)  +  2y8C(a)cosh(w3ax)]acosayda  (3.12b) 
Jo 

2G  ~rxy(x>y)  =  J  [ynE(a)e',wilay/P5  +  y12G(a)e'iw3,ay/p5]asinaxda 

+  f  [2y9A(a)sinh(wiOcx)  +  2y1oC(a)sinh(w3ax)]asinayda  (3.12c) 
Jo 

The  above  expressions  are  for  plane  stress  problems.  It  can  also  be  used  for  plane 
strain  problems  with  the  following  substitutions: 


vyx  -  ‘E^yx’ 


vxy  “  "E^xy* 


EyA  =  Ex  and  EXA  =  Ey. 


3.1.2  Boundary  and  continuity  conditions 


At  the  interface(x1=H1,  x2=-H2) 


®ixxCHi,y)-a2xx("H2.y) 

0<y<oo 

(3.13a) 

'clXy(Pi»y)='c2xy(_H2*y) 

0<y<oo 

(3.13b) 

u1(H1,y)=u2(-H2,y) 

(kycbj  or  b2<y<<» 

(3.14a) 

v1(H1,y)=v2(-H2,y) 

0<y<bj  or  b2<y<°° 

(3.14b) 

cyixx(Hi>y)=P3(y) 

bj<y<b2 

(3.15a) 

xixy(Hi>y)=:P4(y) 

bi<y<b2 

(3.15b) 

At  y  =  0  (cracks  normal  to  the  interface) 

v^Xj.O)^  -HjOc^-aj  and  a1<x1<H1 

(3.16a) 

v2(x2,0)=0  -H2<x2<-a2  and  a2<x2<H2 

(3.16b) 

°iyy(xi’0) — Pi(xi) 

-a^^aj 

(3.17a) 

G2yy(x2>0)='P2(x2) 

-a2<x2<a2 

(3.17b) 

o 

ii 

/ — \ 

o 

£ 

X 

(3.18a) 

S? 

v- 

o 

V-/ 

II 

o 

-H2<X2<H2 

(3.18b) 

At  xi  =  0  or  X2  =  0 

ui(0,y)=0 

0<y<oo 

(3.19a) 
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u2(0,y)=0 

0<y<oo 

(3.19b) 

^lxyC^oO^O 

0<y<oo 

(3.20a) 

o 

ii 

^ — v 

o 

0<y<oo 

(3.20b) 

where  p^Xj)  is  the  normal  traction  on  the  fiber  crack  surfaces  (-a^^  and  y=0), 
p2(x2)  is  the  normal  traction  on  matrix  crack  surfaces  (-a2<x2<a2  and  y=0),  p3(y)  is  the 
normal  traction  on  crack  surfaces  along  the  interface  (x^Hj,  b1<y<b2)  and  p4(y)  is  the 
shear  traction  on  crack  surfaces  along  the  interface  (xj=Hls  b^ycbj). 


Define: 


.  .  3vi(x„0) 

^l(Xl) ~  3xj  5 

**>-*£!* 

a 

^(y)  =^[vi(+Hi,y) 

-  v2(-H2,y)],  <t>4(y)  =  ^[Ui(+Hi,y)  -  u2(-H2,y)] 

(3. 21  abed) 

Then,  from  (3.14ab)  and  (3.16ab),  one  can  write: 

<i>l(Xl)  =  0 

a1<lx1l<H1 

(3.22a) 

<1>2(x2)  =  0 

a2<lx2l<H2 

(3.22b) 

u 

/ — N 

II 

o 

0<y<bj  &  b2<y<~ 

(3.22c) 

o 

II 

^ — N 

V — ✓ 

•e- 

0<y<b1&  &  b2<y<°° 

(3.22d) 

Using  eqn(3.9a)  and  (3.12c),  eqns(3.19a,b)  and  (3.20a, b)  are  satisfied  identically. 

Substituting  eqns(3.18a,b)  into  (3.12c)  we  have: 

E«x)  =  -?G(a),  E*(a)  =  >*G*(a) 

Yu  in 

Then  eqns(3.9a,b)  and  (3.12a,b,c)  become: 


4  r°° 

u(x,y)  =  —  [A(a)sinh(w1ax)  +  C(a)sinh(w3ax)]cosayda 
n  Jo 

+  -  p  E(a)[e-|wllay/Ps  -  ^e-|w3lay/p5]sinaxda 
It  Jo  Yl2 

4  r°° 

v(x,y)  =  —  [p7A(a)cosh(w1ax)  +  (38C(a)cosh(w3ax)]sinayda 
n  Jo 


(3.9a*) 
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-  —  f  E(a)[J39sign(w1)e'lwl*ocy/P5  -  ^Upiosign(w3)e''w3*ay/P5]cosaxda  (3.9b*) 

7t  Jo  Yl2 

— — gy--^-gx(x,y)  =  [  ECaJ^e’^i^y/Ps  -  y^e'^s^y/PsJacosaxda 
4Ex  Jo  ri2 

+  f  ^AC^coshCw^x)  +  2Y4C(a)cosh(w3ax)]acosayda  (3.12a*) 

Jo 

—  ^xyvyx).gv(x)y)  _  f  E(a)[Y5e'^wilay/P5  -  yg-^V 1  w3 1 ay/Ps] acosaxda 

Ziiy  j  Jo  Y12 

+  f  [2Y7A(a)cosh(w!ax)  +  2Y8C(a)cosh(w3ax)]acosayda  (3.12b*) 

Jo 

2G  ~  Cxy(x,y)  =  J  E(a)Yi i [e“*wi ^y/P5  -  e'^w3'ay/p5]asinaxda 

+  f  [2Y9A(a)sinh(w1ax)  +  2Y10C(a)sinh(w3ax)]asinayda  (3.12c*) 

Jo 


Applying  eqns(3.16a,b),  at  y  =  0,  one  obtains: 


v1(x1,0)=  f  E(a)[p9sign(w1)-^Li(310sign(w3)]cosax1da  a^^Hj 

7t  Jo  Y12 

v2(x2,0)=  -z.  f  E*(a)[P9*sign(w1*)-^P10*sign(w3*)]cosax2da  a2<x2<H2 
71  Jo  Yl2 


(3.23a) 


(3.23b) 


3v  (x  0)  2  r°° 

<J>i(xi)=  3XJ’  =  ~  Jo  Yi3E(a)ocsinaXlda, 

3vo(x,,0)  2  r~ 

<>2(x2)=  9x*  -  ~  Jq  Yi3*E*(a)asinax2da 

The  inverse  Fourier  transform  gives: 


i  r 

)_  Yl3a  J  ^ 

a)=Yl3*®  J 

VO 


(^^XjJsinaXidxj 


<j>2(x2)sinax2dx2 


(3.24a) 


(3.24b) 
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Substituting  eqns(3.24a,b)  into  eqn(3.9b)  respectively  and  applying 
eqns(3.17a,b),  one  can  derive: 

fai  Mldt  +  f°°  [2y7A(a)cosh(w1ax1)  +  2Y8C(a)cosh(w3ax1)]ada 

t'xi  Jo 


(3.25a) 


Yu*  T2  ^dt  +  f"  [2Y7*A*(a)cosh(w1*ax2)  +  2Y8*C*(a)cosh(w3*ax2)]ocda 

14  J-a,  t_x2  Jo 


J-&2  1  A2  *io 


(3.25b) 


where  Yi  and  Yi*  are  given  in  Appendix  B. 


Substituting  eqns  (3.24a, b)  into  eqns(3.15a,b),  one  has: 

„-lw,lntv/Sr 


7t(1-Vxyvyx)  (  )  =  ^  f'^sinatdttYie-K'ay/Ps  -  Y2^e-lw3l<*y/p5]cosax1da 
2EX  i  x.-m.Jo  Jo  Y13  <12 

+lim  r[2Y3A(a)cosh(w1ax1)  +  2Y4C(a)cosh(w3ax1)]acosayda  bj<y<b2  (3.26a) 

— 2— p,(y)  =  lim  f  r’^sinatdtYnte-lwitay/Ps  -  e-lWjtay/Pslsinax.da 

2GXy  *,->».  Jo  Jo  Yl3 

+lim  f  [2Y9A(a)sinh(w1ax1)  +  2Y10C(a)sinh(w3ax1)]asinayda  bj<y<b2  (3.26b) 
Jo 

Applying  eqns(A.5a,b),  eqns(3.26a,b)  can  be  further  reduced  to: 


+  Y2(?1)riwJv^1 


J,™^)2  Yl2^+(Hrt)2 

+lim  f°°  [2Y3A(a)cosh(w1ax1)  +  2Y4C(a)cosh(w3ax1)]acosayda 
JO 


ft(l-VxyVyx) 

2EY  1 


bi<y<b2 


(3.26a*) 


and 
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f 

a, 


"YiiM) 


2Yi3  Wily)2 


lw,ly 


lw4ly 


(lw,ly)2 


Idt 


^2  +(Hrt)2  ^+(Hrt)2 

+  lim  r  [2y9A(a)sinh(w1ax1)+  2Y10C(a)sinh(w3ax1)]asinayda 

JO 


n 


2G 


xy 


■p4(y) 


bj<y<b2 


Applying  eqns(3.14a,b),  (3.13a,b)  and  (3.21a,b,c,d),  one  obtains: 


ACoOsinhCwjOHj)  +  C(a)sinh(w3cxH1) 

+  A*(a)sinh(w1*aH2)  +  C*(a)sinh(w3*aH2)  =  Rj(a) 
^ACa^oshCw^!)  +  pgC(a)cosh(w3aH1) 

-  p7*As‘c(a)cosh(w1*aH2)  -  pg*C*(a)cosh(w3*aH2)  =  R2(a) 
^ACa^oshCwjaHj)  +  Y4C(a)cosh(w3aH1) 

-  A,1Y3*A*(a)cosh(w1*aH2)  -  >.1Y4*C*(a)cosh(w3*aH2)]a  =  R3(a) 


[YgAC^sinhCwjOtHj)  +  Yi0C(a)sinh(w3OcH1) 

+  +  X2Y10*C*(a)sinh(w3*aH2)]a  =  R4(a) 

where  R^a),  R2(a),  R3(a)  and  R4(a)  are  given  in  Appendix  A. 

Solving  eqns(3.27abcd)  for  A(a),C(a),  A* (a)  and  C*(a),  we  have: 


(3.26b*) 


(3.27a) 

(3.27b) 

(3.27c) 

(3.27d) 


1  R,(a)  .  .  R,(a).  .  .  R,(a)  .  .  R4(a) 

A<“>  =  cosh(ffllaH,)t  -^-g,(a)  +  +  f(a,  “><“> +  f(o)  n>^ 

1  R,(a)  .  .  R,(a).  ,  .  R,(a)  .  .  R.(a) 

C<“>  -  cosh(co3aH,)[  +  +  f(a)m’<“>  +  f(a)  n>(a)] 

1  R,(a)  .  R,(a).  .  .  R,(a)  .  .  R„(a)  .  . 

A™-eo*«*W  1 13(00  +  m  mjCa)  +  m  n’(a)1 

C*(a)  =  coSh(co3*aH2)[lwS4(«)  +  -^V«>  +  f(a)  «.(«> +  f(a)  n4<°0] 

where  f(a),  gi(a),  hj(a),  mj(a)  and  n^tx)  (i=l, 2,3,4)  are  given  in  Appendix  C. 
Substituting  (3.28a, b,c,d)  into  (3.25a, b)  and  (3.26a* ,b*): 


(3.28a) 

(3.28b) 

(3.28c) 

(3.28d) 


26 


[^~+JtK11(x1,t)]<|>1(t)dt  +  f  Kj 2(x j  ,t)4>2(t)<it 
1 

fb2  fb2 

+  I  K13(x1,t)<)3(t)dt  +  K14(xl5t)cj>4(t)dt 
Jbj  Wb; 


(1-VxyVyx) 

2Yi4Ey 


-ajOc^ai 


r 


K21  (x2,t)4>j  (t)dt  + 


±fV 

J-a, 


■  +  7iK22(x2,t)]  <f>2(t)dt 


r 

J  b. 


u2  r°2 

+  |  K23(x2,t)<j)3(t)dt  +  K24(x2,t)(j)4(t)dt 

/bj  •'bj 

(1— V*xvV*vX) 


-a2<x2<a2 


f 

•'-a, 


r 


KsiCy.MjCOdt  +  I  K32(y,t)<J)2(t)dt 

"ai  — a2 


f 

Jb, 


•b,  *b, 

+  |  K33  *(y ,  t)4>3  (t)dt  +  I  K34’(y,t)<t>4(t)dt 


r 

Jbi 


bj<y<b2 


r 

*,ai 

; 

J  b. 


r 

v-a 


K4i(y,t)4)1(t)dt  +  I  K42(y,t)<j>2(t)dt 


+  I  K43'(y,t)<j>3(t)dt  +  K44’(y,t)cj)4(t)dt 

'b,  \ 


(3.29a) 


(3.29b) 


(3.29c) 
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2G 


xy 


T»4(y) 


bi<y<b2 


(3.29d) 


where 

Kn(xi,t)  ■ 
K12(xj,t) : 
K13(xl5t) 
K14(xlst)  * 
K21(x2,t)  = 
K22(x2,t) : 

K23(x2,t)  = 
K24(xi,t)  = 


Yi4tcJ 

J_ 

Yl4*J 


[k1(x1,tx)e"a(Hrt)iVlw1l  +  k2(x1,a)e'a(H1-t)P5/lw3l  ]da 

[k3(x1,a)e'aCH2‘t)P5*/lwi!f:l  +  k4(x1,a)e'a(^2't)P5:t:^w3:i<^da 
=  Z~Z  L  [Ji(xi>a)  +  J2(xi>a)  Icosatda 

ri  a 7£*0 

J0  [J3(xj,a)  +  J4(xj,a)  ]sinatda  (3.30abcd) 

:  ~  ~Jq  [k5(x2,a)e'aOHi-t)P5/lw1l  +  k6(x2)a)e'a^rt)Ps/*w3l  ]da 

— ~ J  [k7(x2,a)e'a(H2-t)P5*/lw1¥l  +  kg(x2,a)e'a(H2-t)P5*/lw3*l  ]da 
Yl4  TC  0 

^r^|0  [J5(x2>«)  +  J6(x2»a)  ]cosatda 
irl  [J7(x2,a)  +  J8(x2,a) ]sinatda 

a 


Yl4 


(3.31abcd) 


Tyr  .  ,  1  ,  H,-t 

K3i(y>1)  “  oS  ■J’Yinw.k/tf 


/Yiix _ HuL 


+  lim  “f  [koCy.XjjaJe'06^!'1)?^^!!  +  kioCy.Xj.aJe^^r^Ps/^l  ]da 

K*o 

K32(y,t)  =  lim  —  f  [k11(y,x1,a)e'a(H2“t)P5*/|wi*1  +  k12(y,x1,a)e'a(H2't)P5*/|w3*1  ]da 

7CJ0 

„  ..  .  1  /•<*>  rch(w,ax,)T  ,  ch(w3ax,) 

K3,  (y.t)  =  to  -Jo  [ch(WiaHy,(y.a)  +  ch(W3(XH,i/1o(y.a)  Icosatda 


__  1  (•<*>  rch(Wiax,)T  cn^woux,; 

k34  ty-tj  “  to  ~jo  tch(WiOHi/..<y.“)  +  ch(w35aHyu(y.a)  l^atda  (3.32abcd) 
Iw,  ly  lw,ly 

v  ,  v  _  Yii  r _ §5 _  _ Ps _ , 

=  *»*  ™KHr,)> '  ^+(H,  V 

+  lim  —  J  [k13(y,x1,a)e‘a(H1-t)P5/lw1l  +  kl4(y,Xj,a)e"a(Hrt)P5/lw3l  ]da 
1C*0 


chCw^ax,). 
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K42(y,t)  =  lim  -  f  [k15(y,xlJa)e_a(H2't)P5*/|wi*1  +  k16(y,x1,a)e_a(H2_t)P5*/|w3*1  ]da 

w -•»  Ls lcosrada 

**•■>  -  +  ^  (3.33abcd) 

The  derivation  of  the  terms  k[  and  Jj  (i=l  to  16)  are  given  in  Appendix  C. 

Note  that  for  t  =  y  the  integrals  (3.32c, d)  and  (3.33c, d)  are  divergent.  These 
divergent  parts  must  be  studied  and  separated  by  analyzing  the  asymptotic  behavior  of 
the  integrands. 

After  separating  the  singular  parts  in  (3.32c, d)  and  (3.33c, d),  eqns(3.29c,d)  are 
finally  reduced  to: 


f 


P 


K^y.tXMOdt  +  |  K32(y,t)<>2(t)dt  +  y<j)3(y)  +  |  K33(y,t)<J)3(t)dt 

b, 


; 


P2_ 

27C 


f  (M+^(t)d,+  [  KM(WX>4Wdt 

Jbj  Vbj 


(i-vxyyy^3(y) 


2E, 


bi<y<b2 


(3.29c*) 


J  K41(y,t)(j)1(t)dt  +  J*  K42(y,t)<j)2(t)dt+^J  (^+”)4>3(t)dt 


f 


K43(y,t)<))3(t)dt  -  y<|>4(y)  + 


I 


K44(y,t)<|>4(t)dt 


1 


2G 


xy 


■p4(y) 


bl 

bj<y<b2 


(3.29d*) 


where  plf  p2,  p3  and  p4  K33(y,t),  K34(y,t),  K43(y,t)  and  K44(y,t)  are  defined  in 
Appendix  D  The  additional  single-valueness  conditions  for  each  crack  configuration  will 
be  discussed  later. 
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3.1.3  Asymptotic  Analysis  of  the  kernels 

Depending  on  the  crack  geometries,  some  or  all  of  the  kernels  might  be 
unbounded  as  a  — » «*>.  The  unbounded  kernels  will  affect  the  values  of  singularities  at  the 
crack  tips.  Therefore  an  asymptotic  analysis  is  necessary  to  determine  these  singularities. 

Generally  a  kernel  can  be  expressed  as: 

Kij  =  Kijs  +  Kijf  (ij  =  1,2,3, 4)  (3.34) 

where  Kjjs  is  the  unbounded  part  and  Kjjf  is  the  bounded  part. 

Let  k^  (xj,a)  (i=l-8)  be  the  asymptotic  part  of  k^Xj.a),  when  a  ->  then  we 


have: 


Kns&V)  [k,.(x1.a)e^»-*'1  +  k2. ]da 


Yi««J 


=  A-  f"cosh(w1ax1)e‘PttHl~1)Pi;IWl'+H|lw'llda 


vtay»K  J“ 

+  _JAs_  f~cosh(w1ax,)c'at<H,"t*i/lw>h'H,lw'IIda 
^Yl3Yl4^80  ° 

+  tix  ■  f” cosh(w3axi)e'a[CH,",)Pl/IWjl+H'lw,llda 

_  (Ht  —  t)3s/lwi  1+Ht  lwt  I 

7t  85  [(Hi  -  t)P5  /  IWj  I  +Hj  i  Wj  I]2  -  (w^  )2 

K  (Hi  -QP^/lw,  l+Ht  lw3 1 

K  86  [(H!  -t)p5/IWj  l+Hj  lw3l]2  -(w3X!)2 

1  (H1-t)Ps/lw3l+H1lwil 

71  87  [(Hj  -t)p5/lw3  l+Hi  lw1G2“(w1x1)2 

1  (Hi-Ops/lwal+HJwgl 

jc  88  [(Hj  -  t)p5  /  lw3 1  +Hj  I  w3 1]2  -  (w3Xj  )2 


(3.35) 
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Similarly  one  can  derive: 

„  (  u  (H2-t)p;/iw;i+H1iw1i 

K12S(x1#t)  ^  [(h2 _t)p*/|w*I+H1lw1l]2 -(WlXl)2 

v  ii  (H2-t)fc/lw1>l+H1lw3l 

%  t(H2  -t)p;/iw;i+H1iw3i]2  -(w3Xl)2 

1  (Ha-QPl/lwll+H^wJ 

K^5  [(H2  —  t)P; /I wl I  +Hj I Wj I]2  —  (WjXj )2 

1.  (H2-t)p;/lw3,l-t-H1lw3l 

%^6  [(H2-t)p;/iw;i+H1iw3i]2-(w3x1)2 

1  (Hj -t)P5/IWjl+H2lw|l 

K!1S(x2,t)  =  -J,„1[(Hj_t)Pj/|Wi,+Hi|w.|J1_(w.Xi)1 

,  K  (Ht  — t)Ps/lw,l+H2lw3l 

it  ^[(Hj-Opj/lWjl+HalwJIf-CwJxj)1 

1.  (Ht  -t)P;;/lw3l+H2lw*l 

7C  103[(H1-t)p5/lw3l+H2lw;i]2-(w;x2)2 

1  (H1-t)P5/lw3l+H2lw;i 

n  104[(H1-t)P5/lw3l+H2iw;i]2-(w;x2)2 

1,  (H2-t)P;/lw*l+H2lw*l 

K22S(x2,t)  -  s  ^i«[(h2  -t)p;/iw;i+H2iw;if  -  (w,-x2)2 

j_  (H2-t)p;/iw;i+H2iw;i 

jc  no[(H2-t)p;/iw1‘i+H2iw;i]2-(w;x2)2 

1  (H2  -t)P;/lw3l+H2lw*l 

lu[(H2-t)p;/iw;i+H2iw;i]2-(w;x2)2 

i  (H2-t)p;/iw;i+H2iw;i 

7i  112[(H2-t)p;/iw;i+H2iw*3i]2-(w;x2)2 

The  asymptotic  expressions  for  K31(y,t),  K32(y,t),  K41(y,t)  and 


(3.36) 


(3.37) 


(3.37) 
:42(y,t)  are 


obtained  as  follow: 
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HrL 


/Illv 


HrL 


KmM = w + 

P5  P5 

+  Bm  f  f"  D=9-  (y,x1,a)e■“^H'-'>f•"',  +  k„.  (y,x„a)e' ]da 

x,— >h,  7£"v 

-] 


1 


2Ypii['TldWily)2jfTI  tl, +  T2(Y,2)(|w1lyj2 

P5 


P^+Oi.-O2  '““^h-ovo* 


+  Bm  rcos(ay)e-“KH'-0S-,l*'™'"'),-"da 

^Yl3^80  ° 

lim  rCos(ay)e-“I"I'-0S‘'l”'*<H'',')l*‘,da 

wyuX,-  Jo 


— 1 — T  y - 5li - +  y  p-H\ - - ] 

^  ^+(H,-t)2  ^+(H,-,)2 


HirL 


+  P5 


P52 


(H^QPs/lwJ 

12  ,  .,2 


+  P6 


P52 


ff^-Qp  s/lw3l 

l2  ,  .,2 


[(Hj-tjPj/IWjIf+y2  rt>  [(Hj-t)p5/lw3l]2  +  y 


K32S(y.t)  -  P7 


(H2-t)p;/lw;i 


[(H2 -op'/lwfl)2  +  y2  +Pl  [(H2 - t)K/lw;i]2  +  y 


(H,-t)p;/iw;i 


IWiiX 

K41s(y,t)  =  2y|’ti[(1w,Iv)^5~ 


lEilZ 

_A_ 


P9 


PlOT 


[(H1-t)ps/lw1l]2+y2  rl“[(H,  -t)ps/lw3l]2  +y2 


K42S(y.»  'p"[(H2-t)p;/iw;i]2+y2  "Pl2[(H2-t)p;/iw;i]2+y2 
Let  Jj^  (xj,a)  (i=  1-8)  be  the  asymptotic  parts  of  J^x^a),  as  a  -)«, 

have: 

Ki3S(xi,t)  =  -~JT  [Jlm  (xlfa)  +  J2„  (xj,a)  ]cosatda 
T14710 

=  Y 7(^73  +,^74-tM  cos(at)eawi(H,'x,)da 

WuK  Jo 


(3.38) 

(3.39) 

(3.40) 

(3.41) 
then  we 
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+  +^-62  +^g).  f°°  cos(at)e'“W3(Hl  Xl)da 

TTY 14  ^-80  ° 


w,  (H,  -Xj)  _  w3(H1-x1) 

—  Pl3~  ~  \t2  .  .2  +  P 14 


2  .  *2 


‘KCHj-x^f  +  t2  [WaCHj-Xj)]  +t 

Using  the  same  technique  we  can  derive: 

t  ,  i 

K14s(Xi,0  -  Pi5 K(Hi  _  Xi)12  + 12  P>« [Wj(H,  -  X,)]2  + t2 

w* (H2 -x2)  ,  „  w;(H2-x2) 

K23S(x2,t)  -  Pl7  [w.  ^  )]2  +  ta  P  » [w;  (H2  -  x2)]2  + 12 


(3.42) 


(3.43) 

(3.44) 

(3.45) 


K24S(x2,t)  "  P»  [w-  (Hz  _  Xz + 12  +  p20 [W; (H2  -  x2)]2  + 12 

The  coefficients  Xj  (i=85-112)  and  pj  (j=5-20)  are  given  in  Appendix  B. 

The  asymptotic  analysis  for  kernel  K33'(y,t)»  K34'(y»t),  K43(y,t)  and  K44(y,t)  is  given 
Appendix  D. 

3.1.4  Normalization  of  the  singular  integral  equations 

Using  the  following  transformation  formulae: 

x^r,  t=ajs,  x^r,  t=a2S.  ^  -l<r<l,  -1<s<1 

y = ^(bj-b^r  +  |(b2+bj),  t = ^-b^s  +  J^+bj). 

then  t+y=^(b2-b1)[s+r+b0),  b0=~^^  -l<r<l,  -1<s<1 

Substituting  above  into  the  four  singular  integral  equations,  we  obtain: 

I  f1  [_L_+  (r,  s)]<j>°  (s)ds  +  f  a2k^  (r,s)4>2  (s)ds 

7CJ_1  s-r 

d-vAo 


in 


b2-bj 


£iki(r,S)«(s)ds+^J>:4(r.s)«(s)ds  =  -^^I.1"(r) 


(3.46a) 


f  a,k^,  (r,  s)4»r  (s)ds + —  f  [ - +  tta2k^  (r,  s)]^  (s)ds 

J-i  1  z  7c s  — r 
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+  -W  ka (r, s)<t>| (s)ds + J1 1  ki(r,s)«(s)ds  =  -S-^^p;(r)  (3.46b) 

£  a,k;,  (r,  s)]<t>“(s)ds + £  a2k3°2  (r,  s)JE  (s)ds + &-<(.» + £  kj  (r.s)<b“  (s)ds 

-  Si-  (*  (-J-+ - 5— »J(s)ds+i^>-f  ki(r,s)«(s)ds  =  ^^-p3“(r)  (3.46c) 

2%J-i  s-r  s+r  +  bD  2  *-»  2bx 


|\1k;i(r.s)]«(s)ds+J‘i  3,^(1, s)]«(s)ds+^j‘i(^+^^-*(s)cls 

^^f>«(r,s)<l>;(s)<is-^<|>;(r)+^!'|'k;(r,s)«(s)ds=-^.;(r)  (3.46d) 

2  J  1  £  ^ 

where  superscript  "o"  is  used  to  denote  the  normalized  quantities. 

For  orthotropic  materials,  the  dominant  parts  of  eqns  (3.46c)  and  (3.46d)  are 

coupled.  In  order  to  solve  the  singular  integral  equations  numerically,  it  is  necessary  to 

decouple  these  dominant  parts  of  eqns  (3.46c)  and  (3.46d).  Using  the  technique  described 


in  [22],  first  let: 


e-4. 


Pl 

and  A  =  ^ 


— — i 
2 


i 

2 


_  a^  a2k32  ^  _  b2  bt  k33  k 

ai^41  a2^42.  ^  _^43  ^-44 


1- V  V 

- -2-2Lp“ 

2EX 


then  eqns  (3.46c)  and  (3.46d)  are  combined  into: 


A4  +  B.  f  4*  +  F  f  _Jds_+  |'c9ds+f‘  K<J>ds  =  P 
v  7ti  •’-1  s-r  jri^s+r  +  bo  -1-1 


(3.47) 
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Multiplying  both  sides  of  eqn  (3.47)  by  A"1  gives: 

a  +  JL  f  A~lB(i)ds.+i-  f1  A  lB'(t>-+  f  A"1C6ds+  f  A_1K<))ds  =  A_1P  (3.48) 

jtiJ-i  s_r  7ti  J-i  s  +  r+b0  *’“1  Jl 


0 

-i^ 

0 

•  P2" 
-1— 

P! 

■  d'=  a_1b'= 

-i* 

Pi 

i* 

0 

0 

.  p4 

. 

.  P4 

- 

then  the  eigenvalues  jij  2  of  D  can  be  determined  from: 

D  -  |il  =  0  (3.49) 

giving:  n» -Ml  =  0-4,1  «  =  ±J^  =  ±T 
P1P4  V  PlP4  ^ 

Now  let  R  be  a  square  matrix  such  that: 

DR  =  RA  (3.50) 


where: 


A  = 


t4  0 

.0  -c-1 


where  r\1  = 


— ,  and  ti2  = 

Pi 


Pi 

P4 


,  then  eqn  (3.48)  is  reduced  to: 

1  f.  Avds +i  f  R- A-'B'  Ryds  +  r‘R-.A-.ceds 

j^J-i  s_r  7tiJ-i  s+r  +  b0 

+J1  R-1A-1K(j)ds  =  R"1  A_1P 


Introducing  :  <j)  =  R\j/,  where  y  = 


Va 

IVS, 


(3.51) 


’ajC^r.s)  a2c2(r,s) 
_a1c3(r,s)  a2c4(r,s) 


where 
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a  (*L+*l  5l) 

Pi  P4  V^2 


a2(: 


32  ,  42 


ilc° 


k-0  ik° 
a,(-— 

P4  Pi 


.  IK 

vV 


pi 

k° 


a2(~— - 


P4 

ik° 


■) 


Tl2 


P4  Pi 


(3.52) 


-1  A  -1VT3  -  ^2 _ bj 


RAKR  = 


k°(r,s)  k°(r,s) 
k3°(r,s)  k°(r,s) 


ba-b, 

2 


k°  k 

^33  *• 

Pi  P4 


k)+i(  IK^ 

>4  V  Til  Pi 


+  (k^  Tjj  k43  ^  ^  .  ’Hi  ^k33  ^  k^  ^ 

'  Tl2  p4  Pi  Tiz  P4  V  ^2  Pi  P4 


( 


Tl2  k?4  k43 )  .  ri2  ^k33  ^  k^  ^  ^k33  k^  ^ 


Tli  P4  P4 


Tli  Pi  P4 


Pi  P4 


Tlz  kM  ^  'Hj  k43 

■Hi  Pi  \n2  p4 


(3.53) 


and 


R-1A-1  = 


nij  an  2 

-im!  -m2j 


Pi 


-l. 


P2 


-l 


P3 


PlP3P4 

1 


PlP2P4  P4 

Finally  eqns  (3.46c)  and  (3.46d)  are  reduced  to: 

C"1  f1  v  ds  ri  _  ds 


(3.54) 


- —  +f  a^ (r,s»)“(s)ds  +  f  a2c2 (r,s>()2(s)ds 

j-i  s_r  jt  J-i  s  +  r  +  b0  j- i  j-i 

+  t>?~b-1-  k° (r, s)\j f° (s)ds  +  ba  bl  £ k2 (r, s)\|/° (s)ds 


1— vrov 


=  m, 


*y  y* 


2E. 


p3(r)  +  im2 


2G 


■P4(r) 


(3.55) 


*y 


¥4  (r)  -  ~r  £j  ¥4  (s) — ~  £  ¥3  (s)  g+r  +  b  +£3  (r,  (s)ds + a2c4  (r,  sX>2  (s)ds 

+  bz  bl  k3  (r,  s)y3  (s)ds  +  K  (r,  s)\j f°4  (s)ds 
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(3.56) 


Note  that  the  functions  \j/°(r)  and  y“(r)  are  related  by  the  following  relations: 


irt(r)-|[«W-ifc(r)] 


V:(D  =  5H^BW+«(r)l-  « 


(3.57a) 


(3.57b) 


If  both  materials  are  isotropic,  it  is  found  that  —  =  1,  and  the  above  relations 

IK 


reduce  to: 


V3°(r)=|[C0-)-iC«] 
V4°(r)=|[-^3(r)  +  <t>:(r)]  =  -iVl(r) 


(3.57a*) 


(3.57b*) 


(3.46a*) 


Now  to  expressing  eqns  (3.46a, b)  also  in  terms  of  Y^(r)  and  Y^r),  we  have: 
■if  [-5—+ ratk“  (r,s)]c|>° (s)ds + f a,k”2  (r, s)(>2(s)ds 

+  ^L|‘,  [T.vJCsJds  +  T ,y;(s)]ds  =  pffr) 

2  1 

f1  a1k“1(r,s)^i(s)ds+-f1  [—  +  7ta2k^  (r,  s)](>°  (s)ds 
j-i  Ti’i-i  s  —  r 

MWds  +  T 2v;(s))ds  =  p°(r) 

2  1  2y14ri 


(3.46b*) 


where 


T'=k"’+#» 


T.  —  L- 0  _  i  Hlk° 

1 1  K13  %  _  *14 

V 


T,  =  k°n  + 1  l^-K  and  Tz  =  k^  -  i 

Til 


;  & K 


Til 


(3.58abcd) 
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3.1.5  Determination  of  singularities  for  different  crack  geometry 

The  stress  singularity  at  the  crack  tip  varies  with  the  crack  geometry, 
a)  Embedded  cracks  ( |ax|<  H1,  |a2|<  H2  and  b2>  0 ) 

For  embedded  cracks,  the  crack  density  functions  are  of  the  following  forms: 


«(t) 
V?(t)  = 


F°  (t) 


m= 


F2°(t) 


a+tra-tf  * 

f3°(0 _  vnt)=— 

a+t^d-t/3’  d+t)“di-t)?4 


d+t^d-t/2 

f4°  (t) 


(3.59abcd) 


where  F° (t)  (i=l,2,3,4)  are  bounded  at  crack  tips,  and  because  of  symmetry  and 
02=^2- 

Substituting  eqns  (3.59a, b,c,d)  into  eqns  (3.29a, b,c,d)  and  using  the  following 
formula  [33]: 

¥  (z)=IpM>dt=Ip - - dt 

“  t-Z  (t-Zx)  °(z2-t)P"(t-z) 

Fn(z,)cotQtap)  Fn(z2)cot(7tpn)  | 

(z2-z1)p"(z-z1)a°  (z2-z1)“"(z2-z)p"  n 

where  Gn(z)  is  bounded  at  the  crack  tips. 

Then  eqns  (3.29a, b,c,d)  give: 

F,°(-l)cot(na,)  _  F°(l)cot(Tia,)  feYqgfrJ  + 

af'ir+lf'  (2)*(l-rr  2yKEy  F1' 


terms 


F3°(r)  Cl  |-F3°(-l)cot(7ta3)  _  F3°(l)cot(7tp3) 

(l+r)"3  (1— r)fc  i  (2)^  (1  +  r)”3  (2)“3(l-r)P3 

=  m!p3(r)  +  im2P4(r)  +  bounded  terms  bj  <  y  <  b2 

F4(r)  C"1  rF4°(-l)cot(7ta4)  F40(l)cotQt[34) 

(l+r)a4(l-r)P4  i  L  (2)P4(l+r)“4  (2)“4(l-r)P4 


(3.60) 


(3.61) 


StSSS^Sd  -  WMVbl  ,  -  ^l^p2(r)  +  bounded  terms  (3.62) 
(2)“3(r+l)a3  (2)a2(l-r)“3  27^ 


(3.63) 
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=  -  imjp3(r)  -  m^Cr)  +  bounded  terms  bj  <  y  <  b2 
From  eqns  (3.61)  and  (3.62),  it  can  be  easily  shown  that: 

1  a  1 
ai  =  Pi  =  -,  «2  =  p2  =  - 

These  are  the  known  results  for  cracks  perpendicular  to  the  interface  [23]. 
Similarly,  from  eqns  (3.63)  and  (3.64)  we  can  derive: 

cotOtcg  =  -  iC  =  -  i  ]&&- ,  cot(7tp3)  =  i£=  i  1^*- 

\  P2P3  V  P2P3 

cot(7ta4)  =  i C=  i  J&&-,  cot(^p4)  =  -  i£  =  -  i  l^4- 

V  P2P3  V  P2P3 

Using  the  formula: 


(3.64) 


(3.65) 


_i.  .  1  ,  z  +  i 

cot  (z)  =— log - 7 

2i  z — 1 


(3.66) 


We  then  have: 


1 


1 


where 


«3  =  -  +  io,  P3  =  “  '  i(B> 

a4  =  —  -  tco,  p4  =  -  +ico, 
1 

“=^log(R) 


(3.67abcd) 


Eqns  (3.67a-d)  are  of  the  same  form  as  those  obtained  for  interface  cracks  between 
two  isotropic  dissimilar  media  [24-29]. 

b)  Matrix  crack  touching  the  interface(a1=H1,b1>0  and  a2<H2) 

In  this  case,  as  Xj  and  t  approach  to  Hl5  the  kernel  k^x^t)  becomes  unbounded. 

Let  k°lt(r,s)  be  the  normalized  kernel  of  kns(x1,t),  that  is: 

Ht  (l-slPg/lwJ+H^Wjl _ 


k°  (r,s)=^ 


■85 


[H^l-S^/Iwjl+Hjwjf-H^w^)2 
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Rg-sfo/IWtl+HjlWsI  _ 

+  5L86[H1(l-s)p5/lw1l+H1lw3l]2-H?(w3r)2 

H,(l-s)g</lwal+H1lw1l 

+  ?l87[H1(l-s)p5/lw3l+H1!w1l]2 -H?  (Wlr)2 

Hi  g-s)3s/lw3l+H1lw3l 
+  ^88  [Hl  (1  -  s)P5  /I  w3 1  +Hj  i  w3 1  ]2  -  H2  (w3r)2 


2p5 


[- 


H,(l-s+(l+r)^)  H,(l-s+(l-r)^-) 

Ps  rs 


lwtl 


2ps 


1 

[ - 

Hj(l-s+(l+r) 


IwjIIwA 

Ps 


1 

+ 

HjCl-s+g-r) 


lw1llw3l>^ 

Ps 


^>3' 


"I 


P 


XJw3l 


[- 


2Ps  HiCl-s  +  g+r)—— )  HjCl-s+g  r)  3  ) 

Ps  Hs 


Iw.r 


-] 


(3.68) 


Substituting  eqn(3.68)  into  e,n(3.29a)  and  moving  the  bounded  terms  to  the  right  hand 
side  of  the  equation,  we  have: 


If  {J_+xHl^^t - - — - ~“iw,iV 

xJ-i  s-r  2PS  H,(l-s+(l+r)-f— )  H,(l-s+(l-r)  ) 


XgsIWi1, 


Ps 


+1tH^4 


SpTlHl(l-s+a+r)!^)  '  H,(l-s+(l-r)^) 


-] 


P 


40 


1  2R  ltt JWillWjL  „  ...  ,  JwjIIwjL 

zPs  Hj(l-s  +  (l  +  r) — ^ — — )  HjCl-s+Cl-r) — — ) 

Ps  Ps 


] 


+nH)^![ 


2Ps  LH,(l-s+(l+r)!^-)  H,(l-s+(l-r)^-) 

P5  P5 


-]}(j)“(s)ds 


1  — V  V 

- SL- ^-Pi  (r)  +  bounded  terms 


^Yw-Ey 

Applying  eqns  (3.59a)  and  (3.60)  to  (3.69),  and  setting  app^,  we  have: 


F°(-l)cot(jry)  _  rf(l)cot(7cy)  Xw  lw,l  F°(l)  1  [  1 

(2)Y(l+r)Y  (2)Y(l-r)Y  2  p5  _  .  ,  Jw?T  (1+r)Y  a-r)r 


(2)Ysin(7ry) 


LHs  j 


,  K  !wj 


2  p5 


(2)Y  sin(7ry) 


nT  tl+r)T  (1-DtJ 


IWJIW3I 

"ITj 


^87  IW3I 


F°(l) 


2  ps 


(2)Ysin(r/) 


IWJIW3I 

ITu 


1  1 

' O  +  r)’  (l-r)’J 


X«lw,l  fffl)  r  1  I  1  1 
2  P  r...2T>f  Ln4-r^  n-rPJ 


(2)Ysin(7ty) 


lPs 


(l  +  r)Y  (1— r)Y 


1— v  v 

- Pi  (r)  +  bounded  terms 


^Yw^y 

Multiplying  both  sides  of  eqn  (3.70)  by  (l+r)Y  and  letting  r  — >  —1,  we  have: 


(3.69) 


(3.70) 
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1  Wjl 

1 

_  4.^ 

Iwj! 

1 

Ps 

[wfl 

LPsJ 

y  +A,86 

Ps 

Iw^lWjl 

L  ps  J 

+X 


87 


lw3l 

1 

_ L  1 

lw3l 

1 

Ps 

lw,ilw3l 

L  ps  J 

y  +  ^88 

Ps 

[w321 

LPsJ 

(3.71) 


■  =  0 


This  result  is  the  same  as  that  given  in  [23]. 
c)  H-shaped  cracks  (aj=Hi,  b]=0  and  a2>0  ) 

As  aj-^  Hi  and  by-*  0,  the  transverse  crack  intersects  with  the  interface  crack  and  the 
kernels  Kn,  K13,  K14,  K31  and  K41  become  unbounded.  Let  a  be  the  index  of  singularity 
at  the  point  of  intersection,  then  we  have: 

**>-oSr'  “d 

Substituting  above  expressions  into  eqns  (3.29a, c,d),  we  have: 


1  f”‘  r  1  .+ 


-f 

7C  J-H,  t~Xx 


1 


2Ps  "(H1-t+(H,+x,)^£)  (H,-t+(H1-x,)^) 

Ps  Ps 


] 


+lt^4  1 


2fc  l(H1-t+(H1+x,)l2fe!)  (H1-t+(H,-x1)^Y:a;) 

Ps  Ps 


IWillw3l. 


■] 


+  7C 


^87 1  ^3 1 


2Pj  [(H1-t+(H1  +  x1)^S^)  +  (H1-t+(H1-x,)^^)1 

P5  P5 


LJW3L 

+  7E  ■*  3  [ 


2Ps  ‘(H1-t+(H,+x1)^)  (H,-t+(H,-x,)^)  (H‘  ‘2)° 

P5  P5 
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r  jwJCHrJEil  iwanni-A^  ,  *-3w 

Jo  Lpl3 IWjP^Hj— Xj)2+t2  Pl4 lw3l2(H1-x1)2+t2jt“(b2-t)Pl 


IWaKHi-xJ  n  F3(t) 


ff*r  _ It _  _ i_ _ n__£4iiL_ 

Jo  Lpl5  Iw/  ( H !  - xt)2  + 12  Pl6  lw3l2  (Ht  - Xj)2  + 12  t“(b2  -  t)p< 

_  1  Vxy^xpj(r)  +  bounded  terms 
^714^ 


F4(t) 


dt 


t 


«,  ,  1  . 


Hi-t 


Mu  H,-t 


L„,  Vu^+(Hi_t)2  ♦  Y12 

Ps  P5 


(H^OPs  (H,  -  t)ps 

_ IWjJ _ lw3l  Fi(t)  , 

+  Ps(H1-t)2p2  ,..2  p6  (H,-t)2p2  ,  „2}  (H?-t2)g 


I  Wjl 


lw3l' 


,2  +y 


1  — V  V 

- i5L- ^-p3(y)  +  bounded  terms 


Iwily 

rH>  (  Vn  Ps 

J-h,  l2y137cLlw1l2”2 


lw3<y 

Ps 


~+(Hj-t)2  +— 7!^— +(Hl- 


P^ 


P^ 


y  y  Fj(t) 

Iw/  y  lw3l2  y 

*T''»>*fer<-A*s??5V 


2G 


-p4(y)  +  boimded  terms 


*y 


(3.72) 


(3.73) 


(3.74) 
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Multiplying  both  sides  of  eqns  (3.72-7 4)  by  ta  and  letting  t  ->0,  we  then  have  the 
following  three  algebraic  equations  for  Fj(Hj),  F3(0)  and  F4(0)  with  a  yet  to  be 
determined. 


cot(7ta)  ^  IwJ 

*■  /our  \“  OR 


,  K  K1 


(2Hj)°  2  (35  (2H1)asin(7ta)[^-]“  2  &  (2H1)asin(m)[^|^]a 

P  5  P5 


K  iw3i 


Xoo  IwJ 


2  &  (2H1)“sin(7ta)[^^]a  2  & 

Ps  Ps 


)Fi(H,) 


+  2bPcos(^)[|w1|“  +  lw3|a3F3(0)+  2b|sin(^)[|wil“  +  lw3lJF4(0)  ° 


(3.75a) 


(:V1+ 

2Yi3  IwJ  Yn  Iw3I 


’'tP’WV +p«JF!(ri,°11 - L-isrF-(H>) + ^F=(0) 

Ps  IwJ  P5  lw3l  (2Hj)asin(— b" 


+  -2|-[cot(7ta) +-7-^ — r]F4(0)  =  0 
b;  sin(jta) 


(3.75b) 


+(t^t)“] + ’'[p»irl!<ri)“ + p»  x^(rLi)°1 1 - L-*orF><H*) 

2y13  IwJ  IwJ  P5  IwJ  ps  lw3l  (2H,)acos(— ) 


+  £|-[-oot(xa)+— !— -)F3(0)  +  £iF4(0)  =  0 
bj  sm(jta)  b£ 


(3.75c) 


Since  FjfHJ,  F3(0)  and  F4(0)  are  generally  assumed  nonzero,  therefore  the 


characteristic  equation  associated  with  eqns  (3.75abc)  has  to  be  zero.  Numerical 
calculations  indicate  that  a  is  0  as  expected. 
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3.1.6  Numerical  solutions 

In  this  section,  numerical  procedures  are  outlined  for  solving  the  singular  integral 
equations  for  different  crack  configurations  as  discussed  in  the  previous  section.  The  goal 
is  to  determine  the  stress  intensity  factors  and  the  strain  energy  release  rates  for  each  crack 
configuration,  at  the  crack  tips. 

3.1.6.1  Embedded  cracks 

Figure  3.1  shows  the  embedded  crack  geometry.  In  this  case,  the  normalized 
singular  integral  equations  take  the  following  form  as  derived  earlier  in  section  3.1.4: 

—  f  [—  +  7ta1k°(r,s)](l)°(s)ds+  \\  a2k°2(r,s)<j)2(s)ds 

7tJ-i  s-r  J_1 

+  [Ti¥a(s)ds  +  Ti^(s)]ds  =  ^ V^Vg-Pl°(r)  (3.76a) 

2  _1  -^Yu-ky 

f  atk^  (r,  s)<|)°  (s)ds + -  f  [— + Jca2k^  (r,  s)]<S>2  (s)ds 

j-i  jjJ-i  s  —  j- 

+  b2~-^r  [T2\|/°(s)ds  +  T2\jf3(s)]ds  =  (3.76b) 

2  "  "Y 14  ^y 


V3(r)  +  ^rJ11¥30(s)77-“£1V4(s) 
jji  j-i  s— r  7t  J  1 


ds 


s+r  +  b„ 


+  £  k°  (r,  s)\j/3°  (s)ds + £ 


2 

1— v„v 


=  m 


1 


■  2E 


*y 


'rtW-kJ-f  vj(s)— ■ -kl.  |V°(s) 

7C1  J-1  S-r  7t  J-1 


ds 


s  +  r+b„ 


+  £  K  (r»  s)¥a  (s)ds + £ 


1- V  V 


=  -ma 


*y  y» 


1  2E. 


P3(r)-m2-^-p4(r) 


2G 


*y 


■  ■ +£  a  jCj  (r,  s>|>j  (s)ds  +  £  a2c2  (r,  s»)2  (s)ds 
k2(r,s)\|r°(s)ds 

(3.76c) 

+£  ajC3  (r,  s)4>i  (s)ds  +  £  a2c4  (r,  s>j>2(s)ds 
k4(r,s)Y4(s)ds 


(3.76d) 
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All  the  quantities  are  defined  in  section  3.1.4.  Among  the  four  equations  given 
above,  only  three  are  independent  since  xj/g  (r)  and  xj/^Cr)  are  related  through  eqn  (3.57b). 

The  first  three  will  be  used. 

Equations  (3.76a)  and  (3.76b)  are  Fredholm  equations  of  first  kind.  Thus  Lobatto- 
Chebyshev  integration  formula  is  used.  Equation  (3.76c)  is  a  Fredholm  equation  of  the 
second  kind,  and  thus  a  Lobatto-  Jacobi  quadrature  integration  formula  is  used. 


Let  0? (s)  =  “7==T,  ^>=-1^ 


•s/l-s" 


-d  ^®=7lJ^F=(1-s>“(1+s4d'^,<s) 


here  a  =  -  p3,  £  =  -03,  F3°(s)=£  d^^s)  and  d,  =  (dxl  +  i  d^) 

1=0 


Using  the  methods  described  in  [30,37],  we  have: 
n  w  1 

■ r  +  7talK“  (ri,sk)]F1°(sk)  +  a2WkK1°2(ri,sk)F2°(sk)} 


I 

S  7C  Sk-q 


b2  —  bj 


£  £  n1(ri,syj)WIJ4,P®K(s!j)  +  T,(r,,Syj)W^d1Pr,(syp]=^-f^-p:(I;) 

1=0  j=i 

(3.77a) 


1-v  v„ 


n  w  1 

Y  { aiWkK21  (rjjS^Fj0 (sk)  +-^[ - +  7ta2K^(ri,sk)]F2°(sk)} 

"  ^  Sk-Ti 


^2  ^1 


N  n 


£  £  n2(ri’syj)WyJ-d1Pj<ttt>)(syp  +  T2(ri,syj)WyjdiPi<lP)(syy)]  = P2(Ti) 
1=0  i=i 


1=0  j=l 


j,  vrt)  «  A  w,jd,p!a<l(Sa) 

Tti  1=0  J=1  Sjrj  —  Tyi  ^  1=0  j=l  s^+r^  b0 
+  2  [aici(ri>sk)F°  (sk)  +  a2c2(ri,sk)F2  (sk)]Wk 

k=l 


(3.77b) 
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+ 


^2 

2 


S 


2  [T3(ryi,s,j)Wyjd,P<“®(sjj)  +  T,(ryl^)j)W!jd,Pr,(syj)] 


-mi  2g  yxP3(ryi)  +  m22G  p4(ryi) 

(3.77c) 

where  T3  =  k°  (r^s^),  T4  =  -i  j—  k°  (ryi,syj) 

V  "Hi 

and 

w^1  11 

k  2  n-1 

(k=  1,  n) 

w  -  K 

1  n-l 

(k  =  2, ..n-l) 

s  =  cost 11(15  " 11 1 

51  COSl  (n-l)  J 

(k  =  l,...n) 

r7t(2i-l), 
r  _  cosr_ ^ - 1] 

‘  2(n-l) 

(i  =  l,...n-l) 

...  N![r(p + 1))1  r(ct + N + 2)  (  . 

yl~  r<p+N+2)r(a+p+N+3)  v 

(yj=o) 

...  2“*'N![r(a+l)]Jr(P+N+2)  { 

yj  r(a+N+2)r(a+p+N+3)  ^ 

(yj  =  N+l) 

2a4^+1r(a+N+2)r(p  +  N  +  2)  1 

yj_  (N+l)(N+l)!r(a+p  +  N  +  3)  [P^s*)]2 

(yj  =  1,-N) 

the  abscissa  syj  are  the  roots  of  the  following  equation: 

(i-s/)p,k,m(s,)=o 

and  the  collocation  points  ryi  satisfies  the  equation  of  the  N-th  order  Jacobi  polynomial: 

P^-i’-p-1)(rJi)=0 

Wh  and  Pi(ai?)(syj)  (t=0,..N)  are  the  complex  conjugates  of  Wyj 

and  P/0^^)  respectively. 

At  the  interface,  since  a  and  P  are  complex  numbers,  the  weight  Wyj,  the  abscissa 

syj  and  the  collocation  point  ryi  will  also  be  complex.  The  technique  of  using  the  complex 

weight,  abscissa  and  the  collocation  points  are  not  new.  One  can  refer  to  [36-39]. 
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Note  that  generally  dj  (  f  =  0....N)  are  complex  numbers,  and  equation  (3.77c) 
actually  provides  two  sets  of  equations  (for  the  real  and  imaginary  parts).  Together, 
equations  (3.77a-c)  provide  2(n-l)+2(N)  equations.  Since  there  are  2(n)+2(N+l) 
unknowns  be  determined,  the  following  additional  equations  are  required  to  insure  the 
single  valueness  of  the  solution: 

J1  <(>°(s)ds  =  0,  J‘  4>°(s)ds  =  0  and  £  \|f°(s)ds  =  0 

These  in  turn  provide  four  additional  linear  equations: 

t,  Wkf-"(tt)  =  0 

k=l 

X  WkF2°(tk)  =  0 

k=l 


X  X  =  0  (3.78abcd) 

1=0  j=l 

Note  also  that  the  last  expression  provides  two  sets  of  equations 

Once  the  values  of  F°(tk),  F2°(tk)  (k=l,...n)  and  d,  =  d^  +  id^  (l  =  0,...N)  are 

determined,  the  stress  intensity  factors  can  easily  be  calculated  (see  appendix  E  for  details 
of  the  derivation): 

Defining  the  stress  intensity  factors  as  follow: 
for  transverse  crack  in  matrix  and  fiber 

at  xt  =  ai;  k(ai)  =  Lim  ^(Xj-a^  aly(xl50) 

at  x2  =  a^  k(a2)  =  lim  J2(x2  -a2)  c2  (x2,0) 

x2-*a2 

for  interface  cracks: 


J _ 1_ 

P3  2G*y 


k2(bj)  -  i 


P4 

P1P2P3 


■V  V 
»y  y* 


2E. 


kA) 


at  y  =  bi: 
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=  Urn  (b,  -  yr  (y  -  bjp.  H„y)  -  , 


at  y  =  b2: 


1  1  ■W-i,lA' 


P3  2Gxy 


'  P1P2P3  2EX 


Sir  1  1  ^  /u  i  1  p4  1 _ V*yV>» 


=  Lim  (bj  -  y)"3  (y  -  b2r  [— ^lxy^i-y)  ‘ 

y->b2  p3  20^  1 


P1P2P3  2Ex 

After  some  lengthy  calculations  (see  Appendix  E),  we  have 

“«-SSfOT 


^ixx(Hi,y)  ] 


1  1  -W-iF  ^Vgk.W  =  i/i-7” Fj°(-i) 


P3  2G,y 


P1P2P3  2Ex 


(3.79ab) 


(3.79c) 


1  1  -W  -  i  .  p4_  1  =  -,/T?I?(l) 


P,2G„  ■*-"  MPM  2E 
If  the  materials  of  both  layers  are  isotropic,  it  can  be  shown  numerically  that 


(3.79d) 


J _ 1_  _ 

P3  2G^  { 


P4 


1  — v  v 

ny  y* 


P1P2P3  2Ex 


and  eqns  (3.79c, d)  reduce  to 


k2(b!)  -  ik^)  =2Gxyp3  F30  (- 1 ) 


(3.79c*) 

(3.79d*) 


k2(b2)  -  ik!(b2)  =  F3° (1) 

The  strain  energy  release  rate  at  the  interface  crack  tip  can  be  derived  as  follows: 


“  =  -r2£_[E^k±k?  +  -J_±k>] 


(3.80) 


Ay  -^i~C2  2Ex  p2  20^  P3 

Again  if  both  layers  are  made  of  isotropic  materials,  the  above  equation  reduces  to 


AE  2%  1  1  „  2  .  1  21 


Ay  2G*y  Ps 


(k^+k^ 


(3.80*) 


49 


It  should  be  noted  that  when  upper  and  lower  interface  cracks  meet  (that  is  ^=0), 
eqns  (3.76a-d)  should  be  modified  before  proceeding  to  the  numerical  solution.  First  the 
integral  limits  for  interface  cracks  should  be  extended  from  (0,  t^)  to  (-b2,  bj).  Then,  using 
the  following  relations: 

m = «(-t)  «<0 = -  «(-0 

Ki3(xi,  t)  =  Kj3(xi,-t)  Ki4(xi,  t)  =  Ki4(xi,-t)  (i  =  1,2) 

Kj3(y,  t)  =  Kj3(y,-t)  and  Kj4(y,  t)  =  Kj4(y,-t)  (j  =  3,4) 

eqns  (3.76a-d)  can  be  reduced  to: 

—  f1  [—  +  ra^k"  (r,  s)](|>i  (s)ds  +  f '  a2k“2  (r,  s)<t>°  (s)ds 
7t J-1  s-r  J_1 

+  ^-f‘  [T.^Wdst T,v:(s)]dS  =  k^-pf(r)  (3.81a) 

2  J_1  ^7l4^y 

f 1  a^  (r,  s)<J)°  (s)ds + -  f  [— + Jta2k^  (r,  s)]4>°  (s)ds 

J-l  "  kj-i  s  — r 

+^-f,  [T2V;«ds + t, v;(s)]ds = o.sib) 

2  J 1  2y14iiy 

Vs  (r) + ~r  J-l  ^  ®  7~  + 

+  [k°  (r,  s) Ys  (s)  +  k2  (r,  s)\|/4  (s)]ds 

=  mi~  0VJVy"  p3(r)  +  im2-^-p4(r)  (3.81c) 

2Ex  2Gxy 

vo(r)_L-  f1  Y°(s)—  +  f1  a1c3(r,s>i)i(s)ds+  P  a2c4(r,s»>2(s)ds 
jet J-1  s-r  J“1 

+  ^  f  [kj  (r,  s ) Yj  (s)  +  k°  (r ,  s) \|r4  (s)]ds 

Jit 

=  -imIi^^ps(r)-m2-i-p4(r) 

2Ex  2Gxy 


(3.8  Id) 


50 


The  above  numerical  procedure  and  the  single  valuedness  conditions  can  also  be 
applied  to  solve  eqns  (3.81a-c). 

3.1.6.2  Matrix  cracks  touching  the  interface  (a^Hj,  bx>0) 

When  matrix  cracks  touch  the  interface,  the  kernel  F^Cxj.t)  becomes  unbounded. 
And  the  singular  behavior  at  the  crack  tip  changes.  Now,  redefining  <f>°(s)  as: 


*lW  (l-s2)’ 

the  power  of  singularity  y  can  be  determined  from  eqn  (3.71). 

Since  y  is  generally  different  than  0.5,  in  the  numerical  procedures  described  for 
eqns  (3.78a-c),  we  now  use  the  Lobatto- Jacobi  quadrature  for  integrals  containing  <\>° (s) . 
The  rest  of  the  definitions  and  numerical  calculations  in  eqns  (3.78a-c)  remain  the  same. 
The  modified  version  takes  the  following  form: 
o  w‘  1  ,  .  .  ... 

£  {-^[-^+Jia1K-s(r1,sk)]F“(sl)+Wia,Kf1^ri,s1)Ii”(sk)+!>2WlKI";1(r1,Sl)IJ(St)) 
w  K  Sk~Ti 


V-bj  v 


£  £  [T1(r;,syj)WJjd1P1w>(syj)  +  T.(r:,sjj)W„diP^>(syJ)]  = 


1-V  V 


»y  y*  _0 


1=0  j=l 


Pi(n) 


(3.82a) 


n  W  1 

£  { aj  Wk  (rj,  sk  )F“  (sk )  +— [ — —  +  Tta^r^F’^)} 

mmm  «TT’  C  — -  T* 


k=l 


7t  S,-^ 


+  \^-£  £  rrMWriSi)  +  T2(ri,syj)w^d,prf)(syj)]  = 1  ' »  pfo) 


1-V„V„ 


1=0  j=l 


^Yl4^y 


(3.82b) 


rf  f  Wyjd^cs^  f  w^diPrks,) 

h  H  s^-ryi  TC  s  H  S»+r>d  +b„ 

+  S  [alCl (ryi> sk ) Fj° (s’k ) Wk  +  a2c2(ryi,sk)F2°(sk)Wk] 

k=l 
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i  [T3Or^W^If»Cs^+T4(^W1,d.Pr«g] 

^  1=0  j=l 


=  mi ~~  2V|---P3(ryi)  +  to2^-p4(ryi) 


where 


(3.82c) 


k  r(n  +  2-2y)  1  V 

.  21-’-»(n-l)![r(l-7)f  1 

k  r(n  +  2-27)  [PirTl(sktf 

the  abscissa  s'k  are  the  roots  of  the  following  equation: 


(k=l,n) 


(k  =  2,..n-l) 


a-s»2)i'i7J',l(sL)=o 

and  the  collocation  points  r-  are  the  roots  of  the  following  equation 

p£lor~v;)=o 


The  methods  of  evaluating  the  weight,  abscissa  Wk,  tk  for  Lobatto-Chebyshev 
integration  rule,  Wj5  syj  for  Lobatto- Jacobi  integration  rule  and  their  collocation  points  r; 
and  ryi  will  be  the  same  as  those  described  in  section  3. 1.6.1. 

The  single  valuedness  conditions  become: 


I 


k=l 


w;if(tk)=o 


£  WkF°(t*)  =  0 

k=l 

X  X  Wjd^CSyj)  =  0  (3.83abcd) 

1=0  j=l 

The  egressions  for  stress  intensity  factors  and  strain  energy  release  rates  given  in 

the  previous  section  can  still  be  used  in  this  case  except  for  k(aj) 

Redefining  k(ax)  =  Lim  2Y(x2+H2)vc2y(x2,0) 

*-H2 


we  have  (see  Appendix  E): 
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k(a,) 


(H^y^F’CI)  r,  Kl/Ps  .  *  H/Ps 
a- v>;)Sin  w  ^iwj/pJ  102nw;|k|/p5f 


(3.84) 


3.1.6.3  H-shaped  cracks  (a^Hp  bx=0) 

When  at  =  H2  and  b2  =  0,  the  transverse  matrix  cracks  intersect  with  the  interface 
crack  to  form  H-shaped  cracks  as  shown  in  Fig.  3.2.  This  configuration  has  been  observed 
during  the  tensile  test  of  the  uniaxial  ceramic  matrix  composites.  Since  the  power  of 
singularity  at  the  point  of  intersection  is  zero,  which  has  been  confirmed  numerically,  (j)"  (s) 
will  no  longer  be  singular  at  this  crack  tip.  Also,  the  kernels  Kn,  K13,  K14,  K31  and  K41 

become  unbounded.  In  this  case,  because  b2  =  0,  the  numerical  calculations  will  be  based 
on  eqns  (3.81a-c)  instead  of  eqns  (3.76a-c). 

First,  the  singular  parts  of  the  kernels  Kn,  K13,  K14,  K31  and  K41  in  eqns  (3.81a-c) 
will  be  extracted  and  integrated  in  closed  form.  Then  the  modified  equations  take  the 
following  form: 

—  f  +7ca1k1°ls(r,s)](t)°(s)ds  +  f  aj[k“  (r,s)-k°  S(r,s)]<])°  (s)ds 

n;j-i  s~r 

+  £  a2kj,(r,s)(t.;(s)ds  +  i-fi  [T,sV;(s)ds  +  flsv;(s)]ds 

+  [(T.-T^v^s)^  +  (T,  -  T„)^(s)]ds  =  (3.85a) 

2  2y14iiy 

f  a1k“1(r,s)(j)j(s)ds  +  -  f  [— +  7ca2k^2  (r,s)] c}>° s)ds 

J-l  s  — r 

1  *  * 

+  vl',  rr2v:(s)ds  +  T2^(s)]ds  =  -^^-p;(r) 

2  Zy  j4iiy 


(3.85b) 
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\|/°(r)+—  f  \|^(s)—  +  f  a1cls(r,s)(()i  (s)ds  +  f  ajc^r.s)  -  cls(r,s)]<j>°(s)ds 
7ti J_1  s— r  •'~1  1 

+  £  a2c2(r>s)(J)°(s)ds  +  y  £  [k°  (r,  s)^ (s)  +  k° (r, s)V4 (s)]ds 

1 


1— v„,v 


=  m, 


S__2_ 


2E. 


p3(r)  +  im2— -p4(r) 


2G 


*y 


(3.85c) 


where 


and 


T  =  k°  +  i  Hi-  k° 

Ms  K13s  _  *"14s 

V^i 


T  —  V°  i  Hi  tr° 

1  Is  ~  M3s  \  „  M4s 

IK 


:,.(r.S)  =  [Su.  +  iik|i] 

Pi  P4  V% 


Now  defining  (J)°(s)  =  F°(s),  the  definitions  for  the  rest  of  the  displacement  density 
functions  remain  the  same  as  in  section  3. 1.6.1.  Then  we  have: 

«  w’  1 

£  (  +7la,K-„(ri.Sl)]ir(s1;)  +  aIWk[K“(r1,Sl)-  K?u(rllSll)]Iir(Sk) 

S  *  Si-1; 

ajWkK^^F-Wl+^-l;  £  rr,!(r,sypWM4P«)(Syp  +  T1,(r,,s>pW»5,Pr\V] 

*  1=0  j=l 


^2  2  {T1(ri,Syj)-T1((r1,syj)]Wn-dlIf<,Ji)(syj)  +  [TI(ri,skj)-T1I(ri,sXj)]W»diPiaW(syj)} 

*•  1=0  j=l 


1— V  V 


(3.77a) 


.  ,  .  W  1 

£  { at  Wk  (rs  ,sk)Fj° (sk)  +-^[ - +  7ta2K^2(ri,sk)]F2°(sk)} 

S  7C  Sk-Ii 


N  n 


+  *•£  £  [T2(ri,syJ)Wyjd,P®K(syJ)  +  T^pw.a.pr^]  =  % ’.y  ^fri) 


1  •  * 
l-v„v„ 


1=0  j=l 


14 


-2  2 


1=0  j=l  ^yj  *yi 


2  2 


^  1=0  j=l  ®yj  "f  ^yi  "f  b0 


(3.77b) 
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+  2  WJ^a1cls(ri,s,k)F1°(sk)  +  2  W’a^ c^q,^)  -  cls(ri,sk)F°(sk) 

k=l  k=l 


+  Z  WAc2(ri,s1)I?(St)+^X  £  [T,(r(1,syj)W„<i1PlM,(s(J) 

k=l  "  1=0  j=l 


+  T4(r(1,s>j)W!jdiPfS>(sJj)]  =  m1- 


-V  V 

*y  V 

2E. 


p3(ryi)  +  im2^-p4(ryi) 


(3.77c) 


where 


Wk 


2 

n(n-l) 


(k=l,n) 


w-  2 _ 3— 

k  n(n-l)  [Pn_i(sk)]2 


(k  =  2,..n-l) 


the  abscissa  sk  are  the  roots  of  the  following  equation: 

(l-Sk2)Pn-2(s'k)  =  0 

Wk  and  sk  are  the  weights  and  abscissa  of  the  Lobatto-Legendre  integration  rule. 
The  collocation  points  q  are  the  roots  of  the  following  equation 

pn(:l-~1)(r:)  =  o 


Again,  the  methods  of  evaluating  the  weight  and  abscissa  Wk,  tk  for  Lobatto- 
Chebyshev  integration  rule;  Wj,  syj  for  Lobatto- Jacobi  integration  rule  and  their 
collocation  points  q  and  ryi  remain  the  same  as  those  given  in  section  3. 1.6.1. 

The  remaining  four  additional  equations  take  the  following  form: 


Fl°(-1)  =  0 

(3.86a) 

F°(1)  =  0 

(3.86b) 

£  <l>2  (s)  ds  =  0 

(3.86c) 

and 

J_h  <i>i(s)ds  =  j^cMslds 

(3.86d) 

The  first  two  equations  originate  from  the  requirement  that  <J)°(s)  remains  finite  at 
both  crack  tips.  The  third  equation  is  deduced  from  the  single  valuedness  condition  for 
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embedded  cracks  in  the  fiber  if  ^  is  not  zero.  And  the  last  equation  states  the  condition  of 
continuity  of  displacement  at  the  point  of  intersection. 

Eqns  (3.86ab)  can  be  direcdy  used  in  the  numerical  calculations.  Eqn  (3.86c)  can 
be  written  as: 


jr  WkF2°(tk)  =  0  (3.86c*) 

1=1 

Using  the  fact  that 

J^2  <{>3  (s)ds  =  -y-J*  <{>3  (s)ds 

=  -1^-/^  ¥3°(s)ds 
2  °cosh(7tco)’ 


we  have 


.  b,  j  n 

i— dn - 

2  cosh(jtco) 


-H 

Since  |  1  <})1(s)ds  is  real,  therefore, 

J-H, 

Re(do)  =  0  (3.86d*) 

This  result  gives  the  last  equation  that  is  needed  for  the  numerical  calculations. 

The  calculation  of  the  stress  intensity  factors  at  the  embedded  crack  tips  remain  the 
same,  and  needs  no  further  elaboration. 

3.1.7  Results  and  discussion  for  singular  integration  formulation 

It  is  well  known  that  the  interface  bonding  strength  contributes  to  the  toughness  of 
ceramic  matrix  composites.  A  "weak"  bonding  strength  generally  improves  the  toughness 
of  the  composite.  A  direct  indication  of  weak  bonding  between  the  matrix  and  the  fiber 
interface  is  the  existence  of  interface  crack  or  defects  either  due  to  manufacturing  or  due 
to  debonding  during  loading.  In  the  experiment,  it  was  observed  that  under  tensile  loading, 
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matrix  cracks  generated  first  and  then  propagated  to  the  interface.  Once  the  matrix  cracks 
reached  the  interface,  they  did  not  penetrate  through  the  fibers.  This  means  the  fracture 
energy  available  at  the  crack  tip  is  less  than  the  fracture  energy  of  the  fiber,  and  is  larger 
than  the  fracture  energy  along  the  interface  due  the  weak  interface.  Thus,  the  energy  is 
released  along  the  interface  causing  interface  debonding  or  generating  interface  cracks. 
The  actual  role  that  the  interface  plays  has  drawn  a  lot  of  attention  by  both  involved  in 
research  and  manufacturing  of  ceramic  composites.  One  can  refer  to  [45,46]  for  detailed 
discussions.  The  goal  here  is  to  explain  or  simulate  the  damage  mechanism  of  ceramic 
matrix  composites  by  analyzing  the  interactions  between  matrix  and  interface  cracks  using 
different  crack  configurations. 

The  numerical  calculation  for  embedded  crack  configuration  has  been  carried  out 
first.  To  check  the  correctness  of  the  numerical  procedure,  some  simple  crack 
configurations  have  been  studied  first.  For  periodic  transverse  cracks  in  the  matrix  and  the 
fibers  without  interface  cracks,  the  results  match  those  given  in  [23,44];  For  interface 
cracks  without  transverse  cracks,  if  both  Hj  and  H2  are  set  very  large  and  the  materials  of 
both  layers  are  isotropic,  the  problem  is  reduced  to  an  interface  crack  in  an  infinite  domain 
and  the  results  from  [25,43]  for  stress  intensity  factor  at  the  interface  were  recovered. 
Table  3.1  gives  the  results  of  stress  intensity  factors  at  the  interface  for  four  different 
materials  combinations.  As  can  be  seen  from  the  table,  the  two  results  are  almost  identical. 
Also,  the  strain  energy  release  rates  at  the  interface  for  the  above  four  material 
combinations  were  calculated  using  eqn  (3.80*).  For  isotropic  materials,  the  exact 
expression  of  the  strain  energy  release  rate  at  the  interface  was  first  given  in  [40]  as 
follows: 

AE.tt  (n.  +  K.MH.+KfllL-  (k?+k2}  (3.87) 

Ay  2  |11h2[(1+k1)h2  +  (1+k!)(11] 

where  and  } i2  are  the  shear  moduli  of  the  fiber  and  the  matrix  respectively. 
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Table  3.1  Comparison  of  stress  intensity  factors 


Materials 

K,(b) 

K,(b) 

<j0Vb 

o0Vb 

Layer  1 

Layer  2 

Results  from 

Results  from 

[251 

T251 

eqn  (3.79d) 

Aluminum 

1 

1.0000 

-0.1342 

-0.13421 

Steel 

mm 

1 

1.0000 

-0.1443 

-0.14431 

Steel 

Aluminum 

1 

1.0000 

-0.09158 

-0.091576 

Nicalon 

CASH 

1 

1.0000 

-0.071076 

-0.071074 

The  results  using  eqn  (3.80*)  differed  with  that  given  in  eqn  (3.87)  by  a  factor  of 
V2 .  The  strain  energy  release  rates  using  both  eqn  (3.80*)  and  eqn  (3.87)  are  compared  in 
Table  3.2.  The  error  is  within  3  percent.  Considering  the  fact  that  all  the  derivations  are 
based  on  the  orthotropic  materials,  and  the  isotropic  behavior  was  approximated  in  the 
calculation,  the  results  are  quite  good.  Table  3.3  gives  the  Young's  moduli  and  the 
poisson's  ratios  for  the  materials  used  in  the  calculation. 


Table  3.2  Comparison  of  strain  energy  release  rates 


Materials 

Strain  energy  release  rates 

Layer  1 

Layer  2 

2  AE  .  onN 

- in  (3.87) 

7t  Ay 

V2-—  in  (3.80*) 

%  Ay 

Aluminum 

■Ml 

0.397614xl0"5 

0.402399xl0"5 

Steel 

.  1 

0.384227xl0"5 

0.389954X10"5 

Steel 

Aluminum 

0.239708x10"* 

0.239753x10"* 

Nicalon 

CASH 

0.203262x10"* 

0.208256x10“* 
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Table  3.3  Material  constants 


Material 

Young's  modulus 

E(lbin.  -2) 

Poisson's  ratio 

Aluminum 

IxlO7 

0.3 

Steel 

3xl07 

0.3 

Epoxy 

4.5x10s 

0.35 

Nicalon 

2.83xl07 

0.28 

CASH 

13.2xl07 

0.28 

Next,  the  effects  of  the  transverse  and  interface  cracks  on  the  tensile  damage 
behavior  of  the  ceramic  matrix  composite  are  investigated.  Suppose  there  exist  both  small 
transverse  matrix  cracks  and  interface  cracks  or  defects.  First,  let  the  matrix  cracks  be 
allowed  to  propagate  and  the  interface  crack  be  fixed  at  a  small  value  b/H2=0.01.  The 
stress  intensity  factors  at  both  the  matrix  crack  and  the  interface  crack  tips  are  calculated 
and  plotted  against  the  increasing  matrix  crack  length.  Fig.  3.3  shows  the  mode  I  stress 
intensity  factors  at  the  matrix  crack  tips  versus  the  matrix  crack  length  for  both  Vf  =  30% 
and  Vf  =  40%.  The  stress  intensity  factors  increased  first  and  then  decreased  as  the  matrix 
crack  approached  the  interface.  The  decrease  of  the  stress  intensity  factors  as  the  matrix 
crack  approached  the  interface  is  because  of  the  fiber  constraint  effect.  The  mode  I  stress 
intensity  factors  at  the  interface,  however,  remained  unchanged  at  first  and  increased 
rapidly  as  shown  in  Fig.  3.4  when  the  matrix  crack  approached  the  interface.  Fig.  3.5  gives 
the  strain  energy  release  rate  at  the  interface  versus  the  matrix  crack  length  for  two 
different  fiber  volume  fractions.  The  influence  of  fiber  volume  fractions  can  also  be  noted 
here.  Next,  let  the  matrix  crack  length  be  fixed  (a2=0.9H2)  and  the  interface  crack  length 
be  allowed  to  increase,  then  the  stress  intensity  factors  at  the  matrix  crack  tips  and 
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interface  crack  tips  were  calculated  and  plotted  in  Fig.  3.6  and  Fig.  3.7.  respectively.  The 
mode  I  stress  intensity  factors  at  the  transverse  matrix  crack  tips  exhibited  a  steady 
increase  until  the  interface  crack  length  reached  a  certain  value  (b=H2).  Fig.  3.8  shows  the 
strain  energy  release  rate  at  the  interface  versus  the  interface  crack  length.  The  Mowing 
conclusion  can  be  drawn  from  the  above  two  cases:  If  both  matrix  crack  and  interface 
crack  exist  in  the  composite,  the  matrix  crack  will  first  propagate  because  of  the  increase 
of  stress  intensity  factors  at  tile  matrix  crack  tips  and  the  low  fracture  toughness  of  the 
matrix.  Then  the  crack  will  be  arrested  when  it  reaches  the  interface  because  of  the  fiber 
constraint  effect.  On  the  other  hand,  if  the  fracture  toughness  along  the  interface  is  small, 
die  interface  crack  will  grow  due  to  an  increase  in  the  strain  energy  release  rate  at  the 
interface.  Once  the  interface  crack  or  debonding  length  reaches  a  certain  value,  the  fiber 
constraint  effect  is  almost  lost  and  the  stress  intensity  factors  at  the  matrix  crack  tips  start 
to  grow  again  as  can  be  seen  from  Fig.  3.9  where  b/H2  =  0.55.  Eventually  the  matrix  crack 
will  reach  the  interface  boundary.  This  has  been  observed  during  testing.  Another 
important  role  that  tile  interface  crack  plays,  is  to  prevent  the  matrix  cracks  from 
ending  self-similariy  into  the  fibers.  As  can  be  seen  from  the  test  results,  fibers  do  not 
break  at  the  same  locations  as  the  matrix  resulting  in  longer  fiber  pull-out  length.  The 
variations  of  the  stress  intensity  factors  versus  the  increase  of  matrix  crack  length  for  two 
different  values  of  interface  crack  length  are  shown  in  Figure  3.9a  and  Figure  3.9b  for 
fiber  volume  fractions  of  30%  and  40%  respectively. 
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3.2  Finite  Element  analysis 

The  H-shaped  crack  configuration  is  adopted  and  used  to  study  the  effects  of 
interface  cracking  and  its  progression  on  the  tensile  behavior  of  ceramic  matrix 
composites. 

3.2.1  Description  of  the  H-shaped  crack  model 

From  the  experiments,  it  is  observed  that,  after  the  appearance  of  initial  transverse 
matrix  cracks,  the  stress-strain  behavior  starts  to  become  nonlinear.  With  certain  increase 
of  the  tensile  load  (about  10  ksi),  a  regularly  spaced  multiple  transverse  matrix  crack 
pattern  is  formed.  Once  this  stage  is  reached,  further  increase  of  load  generates  few  new 
matrix  cracks  until  final  failure.  Therefore,  during  the  period  from  which  saturated  matrix 
crack  formed  to  eventual  failure,  interface  debonding,  matrix  crack  opening  and  sequential 
fiber  breaking  dominate  the  failure  process.  To  capture  the  effect  of  damage  accumulation 
due  to  interface  debonding,  a  configuration  of  periodic  H-cracks  as  shown  in  Fig.  3.2  is 
used  where  the  interface  crack  length  is  allowed  to  increase  with  increasing  load.  The 
popular  ABAQUS  finite  element  code  is  used  with  strain  energy  release  rate  as  the 
criterion  for  the  interface  crack  propagation.  First,  the  strain  energy  was  calculated  for  a 
certain  interface  crack  length,  then  the  interface  crack  length  is  allowed  to  increase  by 
releasing  the  node  at  the  crack  tip  and  the  strain  energy  is  again  calculated.  Thus  an 
estimation  of  the  strain  energy  release  rate  G  at  the  interface  crack  tip  can  be  determined 
using  the  following  expression: 

O-^  (n=0,l,2, . ) 

where  Un  and  Un+1  are  the  strain  energies  associated  with  the  two  subsequent  crack 
lengths,  AD  is  the  interface  crack  increment  and  t  is  the  thickness  of  the  specimen.  If  this 
G  value  is  greater  than  GCT  (critical  strain  energy  release  rate),  then  the  crack  is  assumed 
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to  propagate.  In  practice,  the  GCT  value  at  the  interface  is  extremely  difficult  to  measure. 
Here  we  estimate  Gcr  by  matching  one  point  of  the  stress- strain  curve  with  the  calculated 
value.  Table  3.4  shows  the  computed  GCT  for  both  Vf  =  30%  and  Vf  =  40%. 


Table  3.4  Estimated  critical  strain  energy  release  rate  Gcr 


Fiber  volume  fraction 

L/2H, 

30% 

2.5 

85.36  (lbs/in) 

40% 

3.0 

62.70  (lbs/in) 

3.2.2  Finite  element  results  and  discussion 

Before  the  finite  element  calculation  is  carried  out,  one  has  to  determine  the 
matrix  crack  spacing  L  and  the  fiber  spacing  2H2used  in  the  model  as  shown  in  Fig.  3.10. 
The  ratio  L/2H2  affects  the  results  significantly.  Matrix  crack  spacing  is  determined  by 
averaging  the  matrix  crack  length  from  the  micrographs  taken  during  the  experiment  It  is 
found  that  the  matrix  crack  spacing  decreases  with  increasing  fiber  volume  fraction  of  the 
specimen.  The  average  matrix  crack  spacing  for  Vf  =  30%  was  found  to  be  L  =  166  pm. 
And  for  Vf  =  40%,  the  value  is  L  =  111  pm.  The  fiber  spacing  varies  with  fiber  volume 
fraction,  the  fiber  diameter  and  how  evenly  the  fiber  filament  is  distributed  in  the  matrix 
media.  The  average  fiber  diameter  for  Nicalon  fiber  is  about  16  pm  as  measured  from  the 
test  results.  Consider  the  case  where  fibers  are  evenly  distributed  in  the  matrix  media, 
then  for  Vf  =  30%,  2H2  =  (0.7)(16)/0.3  =  37.3  pm  and  for  Vf  =  40%,  2H2  =  (0.6)(16)/0.4 
=  24  pm.  Therefore  we  have: 

—  =  1 66/37  =  4.44  for  Vf  =  30% 

2H2 

and  —  =  1 1 1/24  =  4.625  for  Vf  =  40% 

2H, 
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Considering  the  fact  that  fibers  are  randomly  distributed  in  the  matrix  and  that  the 
fiber  diameter  is  not  constant,  the  above  ratios  were  varied  between  2.0  to  3.0  in  the 
calculations.  Fig  3.11  is  a  comparison  of  finite  element  results  with  the  test  results  for  Vf 
=  30%.  Fig.  3.12  shows  the  comparison  for  Vf  =  40  %.  It  can  be  seen  that  finite  element 
results  match  the  test  results  quite  well  at  the  start  of  the  non-linear  stress-strain  curve. 
Then  the  two  results  deviate  somewhat  as  the  load  is  further  increased.  This  is  because  at 
higher  loading  the  fibers  start  to  break  as  can  be  seen  from  the  micrographs  taken  during 
the  tests.  This  contributes  to  the  flattering  of  the  actual  stress-strain  curve.  It  may  also  be 
noted  that  the  finite  element  results  for  Vf  =  30%  match  the  stress-strain  diagram  better 
than  those  for  Vf  =  40%.  This  can  be  explained  as  follows:  with  more  fibers  in  the  40% 
fiber  volume  fraction  specimen,  more  sequential  fiber  breakings  occurred  during  the  final 
stage  of  the  experiment  after  the  saturated  multiple  matrix  crack  pattern  was  formed.  This 
can  not  be  simulated  by  considering  interface  cracks  only.  Fig.  3.13  gives  the  stress 
distribution  in  both  fiber  and  matrix  near  the  interface  crack  tips.  Fig.  3.14  gives  the 
normal  stress  distribution  along  the  interface 
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Fig.  2.2  Specimen  geometry  ( Unit  in  inches) 


Matrix  crack  initiated  Fully  developed  matrix  crack 

away  from  the  defects 

Fig.  2.6  Another  pattern  of  matrix  crack  initiation. 


Fig.  2.7  Micrographs  showing  the  location  of  the  first  matrix 
cracks  initiation  and  their  development. 


Load  =  49.6  lbs 


Load  =  33  lbs 


Fig.  2.9  Further  loading  of  the  specimen  is  possible 

after  failure  for  thin  Nicalon/CAS  II  specimen 


Location  and  the  damage  pattern  just  before 
major  fracture  developed. 


Major  fracture  on  the  specimen 

Fig.  2.10  Brittle  behavior  of  the  thick  Nicalon/CAS  II  specimens, 


Fig.  2.11  Long  fiber  full-out  length  for  Nicalon/CAS  II  with  Vf  =  40% 


Fig.  2.12  Short  fiber  pull-out  length  for  Nicalon/CAS  II  with  Vf=30%. 
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Fig.  2.13  Very  smooth  surfaces  on  the  fibers  being  pull-out, 
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Fig.  2.15  Typical  damage  pattern  for  Nicalon/CAS  II  specimen 
taken  at  different  magnifications:  (a)  300x;  (b)  800x;  (c)  3000x. 
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Fig.  2.16  Typical  fiber  breaking  pattern  for  Nicalon/CAS  II  specimen 


*>ca! 


Fig.  2.18  Matrix  cracks  density  vs.  the  tensile  stress 
for  Nicalon/CAS  II  specimen  with  Vf=30% 


Fig.  2.19  Matrix  cracks  density  vs.  the  tensile  stress  at  room 
temperature  for  Nicalon/CAS  II  specimen  with  Vf=40% 
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Room  temperature 


T  =  250  °C 


T  =  400  °C 


T  =  600  °C 


Fig.  2.26 


Final  failure  surfaces  at  different  temperatures 
for  Nicalon/CAS  H  with  Vf  =  30%. 


(b)  30.3  ksi  (c)  52-°  ksi 

Fig.  2.27  Development  of  matrix  cracks  at  T=600°C 
for  Nicalon/CAS  II  specimen  with  Vf=30%. 


(a)  Pre-existing  crack 
before  loading 


(b)  21.1  ksi 


Fig.  2.29  Development  of  matrix  cracks  at  T=400  °C 
for  Nicalon/CAS  II  specimen  with  Vf=40% 


(c)  56.1  ksi 


(d)  64.1  ksi 


Fig.  2.31  Progression  of  fiber  crack  opening  at  T=400  °C 
for  Nicalon/CAS  II  specimen  with  Vf=40% 


(b)  56.5  ksi  (c)  64.1  ksi 

Fig.  2.32  Progression  of  fiber  crack  opening  at  T=600  °C 
for  Nicalon/CAS  n  specimen  with  Vf=30% 


3.3  Stress  intensity  factors  at  the  matrix  crack 
for  Nicalon/CAS  II  (b/H2=0.001) 


Matrix  crack 


3.4  Stress  intensity  factors  at  the  interface  crack 
for  Nicalon/CAS  II  with  (b/H2=0.001) 


Matrix  crack 
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3.6  Stress  intensity  factors  at  the  matrix 

for  Nicalon/CAS  II  (a2/H2=0.9) 


Interface  crack  length  (b/H2) 


3.7  Stress  intensity  factors  at  the  interface 

for  Nicalon/CAS  II  (a2/H2=0.9) 


Interface  crack  length  (b/H2) 


3.8  Strain  energy  release  rate  at  the  interface 

for  Nicalon/CAS  II  (a2/H2=0.9) 


Interface  crack  length  (b/H2) 


Stress  intensity  factors  at  the  matrix 

for  Niealon/ CAS  II  when  b/H2=0.55 


Matrix  crack 


Stress  intensity  factors  at  the  matrix  crack  tips 
for  Nicalon/ CAS  II  with  Vf=40% 


Matrix  crack 
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Fig.  3.11  Tensile  tests  on  Nicalon/CAS  II 
with  ¥^=30%  at  room  temperature 
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Fig.  3.12  Tensile  Tests  on  Nicalon/CAS  II 
with  V^=  40%  at  room  temperature  (thick  plates) 
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Fig  3=13  Normal  stress  distribution  in  near  field 
with  increasing  interface  crack  length  and  L/2H  =2.5 
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Appendix  A 

From  eqns(1.9ab),  one  obtains: 
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3ui^y>-  =  -  —  f  [ACcOsinhCwiOHj)  +  C(a)sinh(w3aH1)]asinayda 

oy  k  Jo 

+  -f  P  M)[  ^e-iwiiay/ps  +  ^(^e-^s'ay/PsjsinatsinocHjdtda  (A.l) 

71  Jo  Jo  Yl3  Ps  H5  112 

^u2(-H2,y)  _  4_C  (0t)Sinh(w,  +  C(a)sinh(w3*aH2)]asin(xyda 
oy  k  Jo 

■  f  f  -^e'IWl*lay/&*  +  '^*lYl  1  *e-|w3*lay/Ps*]sinatsinaH2dtda  (A.2) 

7tJ0  Jo  Yl3  P5  P5  Yl2 


3Vl^1,y)=  —  f  [ji7A(a)cosh(w1aH1)  +  pgCtoOcosh^oiH^lacosayda 
oy  jtJo 


2f-f 

113  P5 


-  -  r  r 

^Jo  Jo 

+  ■^a^LLpiosign(w3)e“'w3*ay/P5]sinatcosaH1dtda  (A.3) 

Ps  Y12 

^V2("H2,y)=  —  f  [P7*A*(a)cosh(w1*ocH2)  +  P8*C*(a)cosh(w3*aH2)]acosayda 
oy  %  Jo 

_lff 

rcJo  Jo 

+ 1 Q  CTl  1  *Pio*sign(w3*)e-|w3*lay/Ps*]sinatcosocH2dtda  (A.4) 

P5  Yl2 


^%-^rP9*sisn(wi*)e'IWl*l{Xy/*35Si 


Using  the  following  integral  formulas: 


I 


e"atsinbt  dt  =  , , , ,  and  e'atcosbtdt  = 


a2+b2 

equations  (A.1-A.4)  reduce  to: 

au^.y) 
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a2+b2 


a>0 


(A.5ab) 


dy 


=  "-f 

K  Jo 


1  r  ml,  „ 
4  2y^i(' 


[A(a)sinh(w1ocH1)  +  CtoOsinhtw-jaH^asinayda 


(t,y)dt 


(A.l*) 


^U2a^2  =  —  f  [A*(a)sinh(w1*aH2)  +  C*(a)sinh(w3*aH2)]asinayda 

oy  7t  Jo 


u 


(A-2*> 

5Vj-^H— =  -  f  [(^AC^coshCwjaHj)  +  p8C(a)cosh(w3aH1)]acosayda 
oy  %  Jo 

-irsw  (A3iit) 

KJo  ZYl3 


where 


-*rgw  (A-3*) 

^Jo  ZYl3 

1  f  [p7*A*(a)cosh(w1*aH2)  +  p8*C*(a)cosh(w3*aH2)]acosayda 
7t  Jo 

-C^(t-y)d‘  (A-4t) 


r4(t,y)dt 
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—  i  ru  _i_A2 


(A.4*) 


(A.  6a) 


(A.  6b) 


(A.6c) 


Applying  (3.14ab)  and  (3.21cd)  to  eqns(A.l*)-(A.4*): 


4  r 

(j>4(y)  = - [A(a)sinh(w1cxH1)  +  C(a)sinh(w3cxH1) 

7t  Jo 

+  A*(a)sinh(w1*aH2)  +  C*(a)sinh(w3*aH2)]asinayda 

+C 

4  r 

<j)3(y)  =  —  [(37A(a)cosh(w1cxH1)  +  p8C(a)cosh(w3cxH1) 

7U  Jo 

-  p7*A*(a)cosh(w1*cxH2)  -  j3g*C*(a)cosh(w3*aH2)]acosayda 

-if  ^(w)d,+C^°4(,'y)dt 


(A.7a) 


(A.7b) 


Taking  inverse  cosine  and  sine  Fourier  transforms  for  both  sides  of  eqn(A.7a)  and 
(A.7b)  respectively,  one  gets: 


ACoOsinhCv^ocHj)  +  C(a)sinh(w3aH1) 

+  A*(a)sinh(w1*otH2)  +  C*(a)sinh(w3*aH2)  =  Rj(a)  (A.8a) 

and 

P7A(a)cosh(w1otH1)  +  p8C(a)cosh(w3aH1) 

-  p7*A*(a)cosh(w1*otH2)  -  |38*C*(a)cosh(w3*aH2)  =  R2(a)  (A.8b) 

where 

1  f“  2  f“  f*1  d>i(t) 

Rl(a)  =  "  2a  J0  <t)4(y)smaydy  +  —  ■J^G1(t,y)sinaydtdy 

+  raj0  1’  ^G2(W)^ayd«Jy  (A.9a) 

i  f  2  r  r  <f>i(t) 

R2(a)=^Jo  <j>3(y)cosaydy+— -^G3(t,y)cosaydtdy 

“^1  f  ^°4(t'y)cosayd,dy  (A'9b> 


Using  the  following  integrals: 
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J  G3(t,y)cosaydy  =  ^sign(w1)P9(e‘aC^rt)P5/lwi'-e_(X(Hi+t)p5^wil) 

-  rsign(w3)P10^11(e'a(^rt)P5/lw3l  -  e-a(H1+t)p5/lw3l)  (A.  10c) 

2  Yl2 
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eqns(A.9a)  and  (A.9b)  are  finally  reduced  to: 
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Y12 


(A.  lib) 


Applying  the  continuity  conditions  (3.13ab)  and  (3.14ab),  one  derives  the 


following: 
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Jo 
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JnJn  Yl3  r’2 


(A.  12a) 


[°°  2[Y9A(a)sinh(w1ocH1)  +  Yi0C(a)sinh(w3cxH1) 

Jo 

+  X2Y9*A*(a)sinh(w1*aH2)  +  ?L2Y10*C*(a)sinh(w3*aH2)]asinayda 
=  -  f°°  fai^sinatdtYn[e-|wiiay/p5  -  e-|w3 lay/Ps] sinal^ da 

JoJo  Yl3 

-  f°°  h  ^^sinatdtYj ^[eM *  lay/Ps*  -  e-Iw3*!ay/P5*]sinaH2da 

JoJo  Yl3 


(A.  12b) 


Applying  inverve  Fourier  cosine  and  sine  transforms  to  (A.12a,b)  respectively, 


we  have: 

[Y3A(a)cosh(w1aH1)  +  Y4C(a)cosh(w3aH1) 


(A.  13a) 


VI 


-  X1Y3*A*(a)cosh(w1*cxH2)  -  X1y4*C;f!(a)cosh(w3*aH2)]a  =  R3(a) 
and 

[y9A(a)sinh(w1cxH1)  +  y10C(a)sinh(w3aH1) 

+  X2y9*A*(a)sinh(w1*aH2)  +  Vfr0*  C*  (a)sinh  (w3  *  aH2)]  a  =  R4(a) 
where 


(A.  13b) 


2  r  Cn  <bi (t)  2  r  /•_  A,i0o(t) 

R3(a)=-"r  1XTs1(t-y)cosaydtdy+-  00  ai  "^VS2(t,y)cosaydtdy  (A.  14a) 

JoJo  13  -Wo 


2  r  a,  ^i<t)?(0r 
’0^0 


r4(«)=- 

and 


Ml  2 


7c  f°°  fai4y^3(t»y)sinaydtdy ■ 

JoJo  113 


^rS4(t,y)sinaydtdy 


(A.  14b) 


S1(t,y)  =  [-y1 

+  CYi 


H,-t 


Ju, 


Hrt 


Hi+t  .y,,s _ Hi+i 


S2(t,y) '  71 + 


] 


+  [Yi* 


H,+t 


jjwiM 

Ps*2 


v  * 

—h*&d 


H>+t 


■] 


s3(t.y)  =  Yii[(|Wi|V) 


lw1  ly 

Ps 


M 1 

_A_ 


-Tut 


Iwdy 

Ps 


J^IX 

_A_ 


^+(H1+,p_^£+(H1+,)3 


S4(t,y)  -Yh*[(|Wt*|v) 


lwi*ly 
Ps* 


'w»*ly 

Ps* 


Ps*2 


Wo2  lb|S£+(H1-t)2 


-Tn*t 


lwi*ly 

Pl* 


lw,*ly 

Pl* 


(A.  15a) 


(A.  15b) 


(A.  15c) 


(A.15d) 


Vll 


With  the  help  of  the  following  integrals: 


S1(t,y)cosaydy  =  —  ,^j^(e'aCHrt)P5/lw1l  _  g-aCHj+OP^IWjl) 


0 


w3l  _  e-a(H1+t)p5/lw3l) 


S2(t,y)cosaydy  =  —  (e'a(H2-t)(35*/lwj*l  -  e'a@^2+t)Ps*^wi*l) 


0 


+a!lii!&!(e-a(H2-t)|35*/lw3*l  _  e-a(H2+t)P5*/lw3»l) 
2  y12*lw3*l 


f°°S3(t,y)sinaydy  =  ^ —  e_ct(^i+t)P5^wi^) 

Jo  Wl 

-|^i(e-a(Hrt)P5/|w3l  -  e-aCHi+OMws1) 
r°°S4(t,y)sinaydy  =  +2^v'*j  (e-Oc(H2-t)  P5*/l  Wj  *  I  _  e-a(H2+t)P5*/lw1*l) 

Jn  1 


-  e‘a(H2+t)P5*/lw3*l) 


(A.  16a) 


(A.  16b) 


(A.  16c) 


(A.16d) 


R3(a)  and  R4(a)  are  finally  reduced  to: 

R3(a)=^|a‘  <f,(t)[^e-a(Hi-t)P^w1l 

h^f“2  w^.[- to%f«-a(Hj-,)|35*/lw‘*1 +^0fe-a(H2-t)p5*/lw 

13  J-Zz 


and 


13  J-&2 


'a>  ,P5[- j^-je-cxCH.-t'ps/'w!1  +i^p-“CHr')Mw. 

-3. 


(A.  17b) 


viii 


Appendix  B 

Definition  of  material  constants  for  both  layers.  A  superscript  *  will  be  used  for  the 
material  in  second  layer. 

Yi=1  +  2VZi?i  y2=l +  Y3  =  w1  +  vxyp7  Y4  =  w3  +  vyx(38 

P5  P5 

Y5  =  vxy  +  1  Y6  =  vxy  +  3  Y7 =  vxy wi +  P7  Y^  =  VxyW3  +  |38 

Ps  Ps 

Y9  =  -1  +  p7wx  Y10  =  -1  +  P8w3  Yn  =  “  ^  +sign(Wl)p9 

Ps 


Yi2=  r  +sign(w3)pi0 

Ps 

Y13  =  signCwjlPg  -  ^ 
Yl2 

sign(w3)p10 

y  _lrY5._Y6.YlLl 

14  2LYl3  Y13Y12 

5L  K  '-W* 

1  ~  E  l-v*  v* 

x  x  vxy  vyx 

3,-°- 

G* 

^3  =  PsY3  P7Y4 

A4  =  AiY3p7-P7Y3 

A5  =  AiY4PrP.tt 

Jr 

11 

s4 

1 

0 

Ay  =  Yg  “  Y9A2 

Ag  —  Y9 —  Y10A2 

Ag  —  AgAy 

^10  =  ^3^8 

^11  =^3Y9 

A12 = A6  A4 

AJ3  =  A5A5 

^■14  =  Vfe 

^-15  =  ^-6^7 

^16  =  ^3^7 

^17  =  ^3^7  +  ^4p8 

^18  =:  ^4^7 

^19  =  ^3^8  +^-5)^8 

^20  =  ^5p7 

^21  =  ^3  -Y3Ps 

^22  =  Y3?7 

^23  =  P7P3 

^24  ~  ^10^16  “^9^16 

A25  =  ^10^17 

^26  =  ^12^16  “  ^10^18 

A27  =  AgA^  —  A13A16 

A28  =  AgAjg 

A29  =  A12A19  -  A13A17 

A30  =  A13Aj8  —  A12A20 

^31  =  ^lAl6— ^9^16 

^32  =  ^11^17 

A33  =  A16Aj2  -  AuAi8 

^34  =  ^9^22  “  ^14^16 

A35  =  ^-9^21 

^36  =  “Aj4Aj7  -  AjjA^ 

Agy  =  A14A18  ~  A22A12 

Ass  =  ^15^16  —  P7P7^9 

Ajg  =  AgA^3 

^•40  =  ^15^17  “  ^12^23 

^41  =  P7P7^12  —  ^15^18 

^-42  =  ^3^16 

^43  =  ^3^17 

^44  =  ^3^18 

ix 


► 


X-45  =  ^10^16  “  X.nX-16 

X-48  =  ^14^16  —  ^10^2 
X-51  =  ^14^19  +  ^13^21 

X.54  =  -^10^23 

X-57  =  X3X-19 

^60  =  79(^4^19  ~  ^5^17) 

^63  =  —  X-gCYsP  8  ■*"  ^5) 

^-66  =  ^7(P7^19  ~  ^-5^23) 

^-69  =  ^5^17  —  ^4^19 

^71  =  ^a(Pg  (^11  —  ^9)  +  p8^-5(Y9  ~  ^7)) 

^73  =  ^2l(^10  —  X-9) 

^75  =  X-9X.5  —  ^13^3  P  8  (X-gY3  —  X34X3) 

X77  =  P7  (X70P7  “  ^15^3) 

A79  =  Ag(A,4Pg  —  ^5p7  ) 


X-46  —  ^-11^20  ^13^16 

X,49  =  A,13A^2  -  X-^X^q 

^52  =  ^10p7p7  ~  ^15^-16 
X-55  =  ^13^23  “  ^15^19 

X5g  —  —  X,4X,16(X,g  —  y9) 

X-61  =  Y3P7(^10  "  X-9) 
X-64  =  “  P7P7(^10  -  ^9) 


X,g7  —  1  x-^x, 


4^16 


X.47  -  — x,nx,19 

X-so  =  -  X<i0X<2i 

X-53  =  ^15^20  —  X73P7P7 

^56  =  —  ^3^20 

X.59  =  —  X,5X16(y9  —  X^) 

X-62  =  ^(Ya^n  "*■  ^4^2l) 
X-65  =  ^8 (^4^23  ~  P7X37) 


X-68  ~  ^5^16 


^70  =  X-3(P7  (X-10  —  ^ll)  ^  P8^8  —  Y9)) 

X^2 =  P?  (X-nX-5 —  X^X,^)  +  P8  (X^X,12  —  X,jjX,4) 
X74  =  X,i2X<3  —  X.10X,4  +  p7  (X-^Xg  —  X,10y3) 

^76  =  ^23(^10  —  X-9) 

^78  =  Ps  (X35X3  —  P7X9) 

X-80  =  X^4  +  X^5  +  X^6  +  X^7  —  X^g  +  X^9 


X-gi  —  —  X^q  'Kqy  — X^2  +  sign(w1)p9(X73  +  X^4  +  X^5) 


wj  k 


(X-43 —  X,57  +  X^9) 


X-82-' 


X-58  -  X-59  -  X.60  +  sign(w1)p9(X.61  +  X62  +  X,63) 


+  •r5^r(X.64  +  X.65  +  X,66  )  -  ^nP,5  ( X,67  +  X68  +  X,69) 


w, 


w, 


X-g3  =  — (X^o  +  ^7!  +  h,2)  -  ^-sign(w3)p10(  ^3  +  ^4  +  X^) 

Yu  Yu 

_  Ou-Mi^ + v, + X,8)  +  ^4(X„  -  x„ + X,,) 

Yu  kl  Iw3| 

X-84  =  —(^58  +  X-59  +  X-6o)  -  —  sign(w3){310(X.61  +  X,62  +  X,63) 

Yu  Yu 


X 


-  2“-  +  U  +  ^4  <hl  +  >-6!  +  >-69) 

In  M  W 


^-85  —  77^81' 


^■86  ~  78^82' 


^87  —  77^83" 


7l37l4^80  7 13  7 14  A  go  7l37l4^80 

A89  =  -  A70  -  ^71  -  ^72  -  sign(w*)p*  (A73  +  A74  +  ^75) 


^■88  ~  77^84' 


7l37l4^S 


^l7l  |%t(^76  +  A77  +  ^78)  -  ^27ll  %(^43  "  A57  +  A79) 


A^q  —  —  A58  —  A59  A60  sign(w1  )|39  (A6i  +  A62  +  A63) 

_  ^i7*  ot(A64  +  A65  +  A66)  -AjYu  pv(A67  +  A68  +  A69) 

wi 


*  * 

*9!  =  %(A7o  +  A,!  +  Xj2  )  +  %sign(w*)p|0  (A73  +  A74  +  ^75) 

7l2  7l2 

+  A1Y2^1'  %(A76  4-^77  +  Ajig)  +  ^ Y«  t%(A43  “  A57  + 

7 12  w3  |w3| 


713714^80 


*  * 

X92  =  %(A58  +  A59  +  A60)  +  —  sign(w3)P*0  (A61  +  A62  +  A63) 

7i2  7 12 

+  AiY2“^Lp^r(A64  +  A65  +  A66)  +  A2Y11  plr(A67  +  A6g  +  A69) 

7 12  |W3|  |W3| 

A93  =  77^89 -T—  t . -  Ag4  =  YgA^o  *  — T~  A95  =  77A91  .  .  =  y&h)2— 

7 137 14  "”80  7l37l4A,80  ll3ll4/',80 

Ag7  =  -A45  -  A46  -  A47  +  sign(w1)P9(A48  +  A50  +  A51) 

+ ■^[■(^52 + a54 + a55)  -  T11^,5  (A42 + a56 + a57) 

hi  hi 

Ags  =  -  h  -  A32  -  A33  +  sign(w1)p9(A34  +  A35  +  X^g) 

+  T^(A^8  +  ^39  +  X40)  -  T1-1  Y’  (A42  +  A43  +  A44) 


w, 


w, 


xi 


>99  =  “(>45  +  >46  +  >47)  ~ 

Yl2 


— sign(w3)(310(  >48  +  >50  +  >51) 

Y12 


-  In.  +  x54 + XS5)  +  ^(1« + K+K) 

Yu  |w3  |w3 


>100  =  ~ (>31  +  >32  +  >33)  -  ■yLsign(w3)(310(  >34  +  >35  +  >35) 

Yu  Yu 

.  In  Mi(^  +  ^  +  X40)  + 2u%(-  X42  -  X43  +  X«) 

Yu  w3 


W, 


>101  ~  >97  Y7 


YuYm^s 


>103  —  >99  Y7  “7T 

Yl3Yl4>80 

>105  =  —  >45  -  A,46  —  >47  —  sign(Wj  )|39  (>43  +  >50  +  >31) 


>102  —  >98  Y g  *  » 

Yl3Yl4>80 

>104  =  ^lOoYg— M- 
Yl3Yl4^80 


“  >1  Yi  j“£j(>52  +  >54  +  >55)  -^Yu  j~*j(>42  +  >56  +  >57) 
>106  =  —  >31  ~  >32  —  >33  —  Slgn(Wj  ){3  9  (>34  +  >35  +  >35) 

~~  >lYlf^r(>38  "*■  >39  >40  )  ^2  Yll  ^t(~>42  ~  >43  >44) 
Wi  Wj 


*  * 

><107  =  +  >46  +  >47)  +  %sign(w3*)p;0  (X48  +  >50  +  >51) 

Yu  Yu 


+  >lY2^*'  i^*  (>52  +  ><54  +  >55)  +  >zY*i  prr(^42  +  ><56  +  ><57 ) 

Yu  w3  |W,| 


*  * 

><108  =  '^■(>31  >32  >33)  -*■  ^_sign(w3)P10  (>34  +  >35  +  >35) 

Yu  Yu 


+  >1 Y2  1^*  (>38  +  >39  +  >40)  +  >^Yu  (  ><42  >43  ^  >44)  ^ 

Yu  w3 


XU 


^109  ~ ^”105  Y7  *  *  y 

Yl3Yl4^80 

*110  “*106 Y«  *  .  * 

If  13  1 14  A80 

*111  “W/?  •  *  y 

113  J14A80 

*112  “*10878  *  *  y 

Yl37l4A'80 

p  _  73*81  +  74*82 
^13  *80 

p  73*83  +  74*84 

^13  ^-80 

n  __  73*89  +  74*90 

P7  xt-aK, 

n  __  73*91  +74*92 

8  WnK 

P  . ...  79*81  +  Yl0*82 
^13^-80 

n  .  79*83  +  YltA-84 

Pl0_ 

_  79*89  +  7l0*90 

n  __  7 9*91  +  YlO*92 

Pn  vr»K 

Pl2  *a. 

_  'YqO^Ti  +^-74  +^75) 
^14*80 

_  _  Yg(*61  +  *62  +  *ffi) 

14  ^14*80 

P  .  Y?(*70  +  *71  +^72) 

78(*58  +*59  +*6o) 

7ry14X.g0 

Pl6  ^14*80 

_  Y?(*48  +  *50  +  *5l) 

Pl7 

n  _  78 (*34  +  *35  +  *36) 

p,,_  «o» 

_  Y 7  (*45  +  *46  +  *47) 

Pls~  *r»K, 

_  _  7  8  (^"31  +  *32  +  *33) 

P20_  *llK 

Appendix  C 


fj(a)  =  tanh(w*1aH2)tanh(w*3aH2)  f2(a)  =  tanh(w*3otH2)tanh(w1(XH1) 
f3(a)  =  tanh(w*3aH2)tanh(w3OH1)  f4(a)  =  tanh(w*1ocH2)tanh(w3aH1) 
f5(a)  =  tanh(w*1cxH2)tanh(w1cxH1)  f6(a)  =  tanh(w1aHi)tanh(w3otH1) 
f(a)  =  X^f^CX)  +  A^f^CC)  +  A^gf^CC)  +  ^7^ (OC)  -  A^gfj^)  +  ^gfgCCX) 
gi(a)  =  X70tanh(w*3aH2)  +  X71tanh(w*1aH2)  +  X72tanh(w3aH1) 
h^a)  =  X73f1(a)  +  ^74f3(a)  +  ^5f4(a) 
nil  (ex)  =  ^7gfi(cx)  +  A77f3(cx)  +  X^gf4(cx) 

nt(a)  =  A43tanh(w*3aH2)  -  X57tanh(w!i!1aH2)  +  X79tanh(w3OH1) 
g2(a)  =  A58tanh(w*3aH2)  +  X59tanh(w*1ocH2)  +  ?i60tanh(w1cxH1) 
h2(a)  =  X61fj(a)  +  X62f2(a)  +  X63f5(a) 
m2(a)  =  ^64fi(a)  +  X65f2(a)  +  X66f5(a) 

n2(a)  =  A67tanh(w*3aH2)  +  X68tanh(w*1aH2)  +  X69tanh(w1aH1) 

g3(a)  =  A45tanh(w*3aH2)  +  ?i46tanh(w3aH1)  +  X47tanh(w1aH1) 

h3(a)  =  X48f3(a)  +  A5Qf2(a)  +  X51f6(a) 

m3(a)  =  X-s2f3(cx)  +  A,54f2(cx)  +  A,5jfg(cx) 

n3(a)  =  A42tanh(w*3aH2)  +  X56tanh(w3aH1)  +  ?i57tanh(w1aHi) 

g4(a)  =  X31tanh(w*1cxH2)  +  /^tanh^aHj)  +  X33tanh(w3aH1) 

h4(a)  =  ^4f4(a)  +  ?^5f5(a)  +  X36f6(a) 

m4(a)  =  ^38f4(a)  +  X39f5(a)  +  X40f6(a) 

n4(a)  =  -  X,42tanh(w*1aH2)  -  X43tanh(w1aH1)  +  ?t44tanh(w3aH1) 

+  sign(Wj)P9h2(a)  +  ^m2(a)  -  -^sn2(a)  ] 


xiv 


**■•*  ■  tgl(a)  •  ^>(w3)P„Ma) 


Yl1.Y2p5  „  .  Ylllk,  f0c')  l  + 


2cosh(w3aH1)f(a)y13L  Yi2^ 


-  ^sign(w3)P10h2(a)-  ^^-m2(a)  +  -^^2(a)  ] 

w-"»  -JSSwi-w  -  WWW) 


_*VR.* 


*Tt.R  * 


-  sign(w!*)p9*h2(a)  -  !wT|  m2(°0  -  |Wl*|  “2(°0  ] 


+  ^sign(w3*)P10*h2(a)  +Yi^T2|w31*i5_ul2(a)  +Y1|w3*I  li2(a)  ] 

_  y^shfoooQh^a)  Y.coshCw.axQh^a) 

Ji(xi»a)  -  coshCwjOH^fCa)  J2^xi’a;  "  cosh^aH^cx) 


oCoshiWQax, 


T  /„  _  (7— - K?!'"'.  T  ft)  =  -  •“ - ^ - 1 

J3^l’a'  ‘  coshCwiaH^fCa)  15  '  coshCv^aH^fta) 

Ww»  -  iSSj  )WI.,W) 

+sign(w,)P9h4(a)  +-^m4(o)  -  -^^n4(a)  ] 

-^ign(w3)P10h3(a) 

-  ^sign(w3)p10h4(a)-  ^^m4(a)  +  ^n4(a)  ] 


XV 


wwa  ■  J»,SwS„-1',’<*1  ‘  '*“*'**’"'w 


■  R  .* 


*Oi-R  * 


-  sign(w1*)P9*h4(a)  -  J[Lj^pm4(a)  -  |wl|  u4(°0  1 
k^“>  +^n<w,*)P,o*h,(a) 


_*v  *0i  _R_* 


Yl2  lw3  1  IW 

+  ^sign(w3*)Pio*h4(a)  + 
J12 


, -7,-/1,  p<-  yn'/^Ps'  ,, ,  Ys*cosh(w3*ax2) — Yu*  ,, 

V  iw^l  ul3(a)  +  lwi*l  ll3(oc)  1  +  2cosh(w3*aH2)f(a)Y13*[  Yi2*§4 


*v  *0l  .R_* 


*0i  R-* 


fll  f ?.  1r5  ...  (n\  ■  fli  '^3-  (a\  1 

Yl2*  lw3*'  4  lw3*l  U4W  j 


Y7*cosh(w1  *ax,)h3(a)  _  .  .  Y8*cosh(w3*ax?)h4(a) 

J5^x2’a) "  cosh(w!*aH2)f(a)  e(x2>a)  Cosh(w3*aH2)f(a) 


YI*cosh(wl*oex,)g,(a)  YR*cosh(w3*ax2)g4(o 

h\x2’a>  coshCw^ocH^fCa)  8^2,  j“'  cosh(w3*aH2)f(a) 


+^i,n*(w‘ wSn^“n 

+  *4»(a)]  +  gC°t^[NM 

Y12Iw3I  1  lw3l  2f(a)Yi3cosh(w3aH1)  Yi2 

-  ^Sign(w3)p10h2(a)  -  ^^.m2(a)  +  1 


-  * w  -  si8n(w'*)fe‘h'(“) 


XVI 


-SSiSi;' 

-  sign(w,*)p9*h2(tt)  -  -  *  ty^''"2*00 1 


Y.cosa^coshCw^xQ,  111  (a)  +  ^sign(w3*)p1o*h1(a) 
k!2(y»xi»a)  -  2f(a)Yi3*cosh(w1aH1)  L  Y12  Y12 


V  *y Yn*kP<*  >  \  l  _L  YdCosajLcosh(w3otx1)  r3&i!. „  fny 

+  Yi2*Y  ,w3  1  mi(a)  +  ,Wl*'  l  0C  +  2f(a)Yi3*cosh(w3ccH1)  Y12* 

+  ^LLrsign(w3*)p1o*h2(a)  +  '  \  \^  *\a^a^  lw3*l  2 


J5(y'a) =If(^fcosay 

J„(y.«)=-3^fcosay 


J,„(y,a)  =  JiM^cosay 


Ji2(y»a)=" 


f(a) 


f(a) 


cosay 


k13(y,X!,a) 


y«sinay.smh(w,ax,)  t  (a)  +sign(Wl)p^1(a)  + 
2f(a)Yi3  sinh  (w^Hj )  1 


JiA 

IwJ 


n^a)  ] 


Y^may-sinhtw^,)  (o)  +sig„(Wl)p,h2(a) 
2f(a)Yi3  sinh(w3aH1) 


JiA 

Iv^l 


n2(a)  ] 


k14(y»xi,a) 


_  Yosina^ sinh (w^)  III  /  }  _  ^*sign(W3)p10h1(a) 
_2f(a)Y13sinh(w1aH1)  Y12  Y12 


YmSinay  sinh  (w3axt) 
2f(a)Yi3  sinh(w3aH1) 


[*“•&«*> 

Yl2 


-  ^LLsign(w3)|310h2(a) 

Yl2 


IllY2Ps  _  .  lufii 


(a)  +  ■, 

Y12  lw3l  lw3 


n2(a)  ] 


ki5(y»Xi,a) 


Yc)sinay_sinh(w.iax1)  [  (a) .  sign(Wl*)p9*h1(a) 

2f(a)Yi3*  sinh  (w1aH1 ) 


xvu 


k16(y,x 


Ti*xiPs*  (g.T"*^*„  (a) ,  |  JuggLsmh(w3ttx,) 
lwi*l  l(  )  lw,*l  l(-  j  J  +  2ffaYv,,*c,w™,  nU  \ 1  ZlK  } 


2f(cc)Y13*  sinh(w3cxH1) 


sign(Wl*)p,*h2(a)  -  Y'*w1'*|S*m^<<x)  '  r"|w^|5*“2(Ct)  1 


'•K)  ‘  ^g‘(a>  +^ign<W3*>P,0*h‘(O) 


+  ^iSn(W3*)^10*h2(a)  +  y!1*Y2|W3*|5  “2(0c) ] 


.  .  .  Yoh,(a)  . 

Ji3(y>a)  =  f(Jx)”smay 


Ji4(y.cx)  =  Tl^(“"sinay 


Ji5(y>a)=-Y<f(^)X')'sinay 


Ji6(y»a)=-Ylj(2()a)'sinay 


xvm 


Appendix  D 

Consider  the  following  integrals  in  eqns(1.29c,d): 


I 


ibi  rb2 

K33,(y,t)^(t)dt  K34,(y,t)4>4(t)dt 


r 

J  b, 


f 

W  b. 


K43’(y,t)<l>3(t)dt  K44’(y»t)et>4(t)dt 


where 


K 


33 


1  f°°  chCwjaxj)  _  ,  .  ch(w,ax,) . 

'<**>  =  J”  ;  W.aH,)1^ + d^5yJ'“(y'a)  Icosatda 

Jo 


+  chtwjgx,)  liML 
chCwjOtHj)  f(a)  ch(w3ocH1)  f(a)  J  y 


chCw^y^Ca) chCw^ax,)  y4h2(tx)  cosa(t.v)+cosa(t+wma 
chCwjOcHj)  f(a)  +ch(w3aH1)  f(a)  JL  (  y)+cosa(t+y)jda 


K 


.„ir 

7t  I 

Jo 

=  I 

Jo 

J»“ 

ch(w,ax,)  ,  .  ch(w3ax,l  .  . _  . 

tch(w1oH1)I|,fy,a)  +  ch(w3aH1/12  y,a  ]smatda 

0 

ch(w,«x,)r3g,(q)  +  ^(w3ax,)j4g2(a)  Qs  gtda 

7t  I  ch(w1ocH1)  f(a)  ch(w3aH1)  f(a)  J  y 

Jo 

_  _  Um  JL  f  |.ch(w1ax1)y3g1(a)  ch(w3ax1)y4g2(a)  vUdnaft+vUdry 
”  xS,27t  chCwjaHj)  f(a)  ch(w3cxH1)  f(a)  Hsina(t-y)+sma(t+y)]da 

Wo 

1  f  sh(w,ax,)  sh(w3ax,) 

'(y'°  “  J5  ;  [ch(w1cxH1),'2(y'“)  +  ch(w3aH1),»(y'a)  Icosatda 

Jo 


K 


43 


(D.labcd) 


(D.2a) 


(D.2b) 


7 


XIX 


=  umif 

7C  I 

Jo 

=*di  [ 

Jo 


fgMM  sh  w,ax,)r,^(a)  satdcc 

lch(w1aH1)  f(a)  ch(w3aH1)  f(a) 


(mo) 

Lch(w1aH1)  f(a)  ch(w3aH1)  f(a) 


sh(w,ax,)_  .  .  sh(w,ax,) . 

+  ch^aH,)1^  lsma,d“ 

0 


=.toir 

%  I 
*0 

f 

Jo 


sh(w,«x,)K,(g)  y»<%l3ufefia jsfaaysinmda 
^(wjaH!)  f(a)  ch(w3aHj)  f(a) 


[»W|C5^(^w  “5 >/’rr)l[cosa(t-y)-cosa(,+y)]da  (D2.d) 

xhCwjaHj)  f(a)  ch(w3aH!)  f(a) 


As  a  approaches  infinty,  we  have: 


cnCw^x^n^cxj  cniw3uxi;  jAn^  .-aw.CH.-x,)  +  P  e-aw3(HrXl) 
WiOHi)  f(a)  +ch(w3aH1)  f(cc)  J"^e  1  1  u  +  & 


(w,ax, 


ciHw1«x,;T3g1w  cn(w3axi;  r4^u;  .-aw.CH.-x,)  +  p  e-aw3(HrXl) 

W.aH,)  f(a)  +  ch(w3dH1)  f(a)  J"^e  lV  1  U  + 

sh(w,ax|)y,,h,(cO  shCw^x,)?,^^  aw  m.-x.)  +  f  e-aw3(Hrx,) 
‘chCwjaH,)  f(a)  +  ch(w3aHj)  f(a)  1  1  1  +i*e 


JHWjCKlLMiW  Im62W,  ,  e-aw,(H,-x.)  +  f  e-aw,(H,-x.) 

ch(WjOtHj)  f(a)  +  ch(w3aHj)  f(a)  J"?’e  1  1  1  +?»e  3  1  1 

where  £i(i=l-8)  are  material  constants  and  are  expressed  as: 

e  _ _ Y?  (^+^4+^-7';) _ 

1  ^24  +  ^25  -**  ^26  4"^27  —  ^28  +^29 

e  _ _ Y4(Xfii+X,fi7+X.6xj) _ 

2  ^24  +  ^25  +  ^26  +^27  -  ^28  +^29 

e  _ _ Y7  (^70+^7 1+^77.) _ 

^4  +  ^25  +  ^26  +^27  -  ?^g  +^9 

e  Yd^s+^g+^o) _ 

4  +  ^25  +  +X^7  —  ^2g +X^9 

t  _ _ Yq(^77+^74+^7s) _ 

^24  +  ^25  +  K.6  +^27  “  ^28  +^29 


sh(w,ax, 


(D.3a) 


(D.3b) 


(D.3c) 


(D.3d) 


(D.4a) 


(D.4b) 


(D.4c) 


(DM) 


(D.4e) 


XX 


YlP^A?4^) - 

^  +  ^25  +  ^26 +^27  ~  ^28 +^29 

Vq(^n+?t7 1+^72) - 

_  ^24  +  ^25  +  ^26  +^27  “  ^28  +^29 

_ _ Yin(^a+^Q+^6o) - 

_  X24  +  ^25  +  ^26 +^27  ~  ^28 +^29 


(D.4f) 


(D4.g) 

(D4.h) 


Using  the  following  formulas: 


Km  f  e-aw1(H1-x1)[cosa(t-y)  +  cosa(t+y)]da  =  jc[8(t-y)  +  8(t+y)]  (D5.a) 


J*°°  1  l 

e-aw1(H1-x1)[sina(t.y)  +  sina(t+y)]da  =  (t_y  +  t+y) 

0 


The  asymptotic  values  of  eqns(D.2abcd)  become: 


(D.5b) 


Kjn  -L  f  [i  e-aw^Hi-xO  +  ^e-aw3(H1-x1)][cosa(t-y)  +  cosa(t+y)]da 

x,^h,  2ji  I  1 

Jo 


-  —  [5(t-y)  +  8(t+y)] 

2r 


(D.6a) 


-  lim  —  f  £  e'awi(Hrxi^  +  ^e-aws^r^HsinaCt-y)  +  sina(t+y)]da 

x,->h,  2tc  I  3 

Jo 


=  _  p2.(— +  — ) 
2tc  t-y  t+y' 


(D.6b) 


^  _L  f  g  e-aWiCHj-Xj)  +  ^e-aw3(Hrx1)][_sina(t-y)  +  sina(t+y)]da 

-,-H,  27t  I  5 

Jo 


=  £l(_  —  +  — ) 
271 1  t-y  t+y' 


(D.6c) 


-  lim  —  f  [^.e"awi(Hrxi)  +  ^e-^CHi-JcOjfcosaCt-y)  -  cosa(t+y)]da 

2%  I 

Jo 


XXI 


__  p4 


[S(t-y)  -  5(t+y)] 


(D.6d) 


where  Pi  =  +  £2 »  P2  =  £3  + *  P3  =  £5  +  ^6  ^  P4  =  £7  + 

For  0<b!<t<b2 ,  one  has  the  following  relations: 


fb2 


I 


<j)3(t)S(t-y)dt  =  -y  <j)3(y) 


(D.7a) 


Pi 


*{ 

r 


4>3  (t)  5  (t+y  )dt  =  0 


<t>4(t)S(t-y)dt  =  —04(y) 


(D.7b) 


(D.7c) 


^  |  <S>4(t)8(t+y)dt  =  0 


(D.7d) 


By  adding  and  subtracting  the  asymptotic  value  under  the  integral  signs  in 
eqns(D.2a-d),  and  substituting  them  back  into  eqns(D.a-d)  respectively,  we  obtain: 


[ 


K33'(y,t)<|>3(t)dt  =  4>3(y)  +  I  K33(y,t)<|>3(t)dt 

bj  \ 


r 

J  b. 


(D.8a) 


f 

Jbx 

r 

J  b, 


K34'(y,t)<t)4(t)dt  =  -  te¬ 
nt 


;)<t>4(t)dt+  I  K34(y,t)<)>4(t)dt  (D.8b) 


Ko(y,t)««)dt  = 

1% 


Jbj  »b 

[  (-^+^)Wt)dt+  [  K43(y,t)cf>3(t)dt  (D.8c) 

Jb,  »bI 


xxu 


K„Xy,t)<MO*=-y'My)  + 


K44(y,t)<J)4(t)dt 


where 


_l_f  ch(w,  ocx,  >7,11,(00  chCw^ax,)  Yah,(a) 

K33(y,t)-Iim27c  chCwjOcHj)  f(a)  +ch(w3aH1)  f(a) 

Jo 

-  teje-awiCHi-Xi)  +  ^2e-aw3(HrXi)] } [C0Sa(t-y)  +  cosa(t+y)]da 


K3)(y,t)=-limif 

Jo 

-  [^e-awiCHrXi)  +  ^4e-aw3(Hi-x1)]}[sina(t-y)  +  sina(t+y)]da 


ch^axQY^Ca)  chCw^axQYd&CoO, 

ch(wjOcHi)  f(a)  +ch(w3ocH1)  f(a)  J 


Jo 

-  [^e-awiCHi-Xi)  +  ^6e-aw3(Hrx1)] }  [-sina(t-y)  +  sina(t+y)]da 


sh(w1ax1)Yoh1(a)  sh(w3ax,)Ymh?(a) 
ch(Wj(xHj)  f(a)  +ch(w3aH1)  f(a) 


K„(y,t)=-Umij 

Jo 

-  [£7e'awi(Hrxi)  +  ^ge‘aw3(Hrxi)] }  [C0Sa(t-y)  -  cosa(t+y)]da 


.shCw^QYogiCa)  sh(w^ax,)  YinB?(a) 
ch(w1aHj)  f(a)  +ch(w3aH!)  f(a) 


] 


(D.8d) 


(D.9a) 


(D.9b) 


(D.9c) 


(D.9d) 
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Appendix  E  Derivation  of  stress  intensity  factors  and  strain  energy  release  rates 
i)  Case  of  embedded  cracks 

The  stress  intensity  factors  at  the  transverse  crack  tips  are  defined  as: 

at  Xj  =  ai:  k(ax)  =  Iim  a^  cr,  (xx,0)  (E.la) 

at  x2  =  &2:  kCaj)  =  Iim  ^2(x2-a2)  c2y(x2,0)  (E.  lb) 


From  eqns  (3.29a),  for  lxxl  >  aj  we  obtain: 

<VM)  =  — Mldt + a»  (x  o) 

where  a“y(xj,0)  is  a  bounded  function,  and 

tftl  _  m  ms™ 

V?7?  (t-aI)w(t  +  a1)1'2 

From  [33]  if  we  introduce  the  sectionally  holomoiphic  function  cp(z)  as 


then  we  have: 


cp(z)=  — P  ■^^■dt 

K  t  -  Z 


,  x  F,  (-a,  )eilc/2 
9(z)= - - 


*M) 


+  q>o(z) 


(2ai)1/2(z  +  ai)1/2  (2a1)1/2(z-a1)1/2 

Substituting  (E.3)  into  (E.2)  and  from  the  definition  of  (E.la),  we  have: 


k(aj)  =  - 


2r“Ey  1  f ,(»,)- 


Similarly 


k(a2)  =  - 


2y;4e; 


*a-v;v;)^  ~  n(i-v;v;) 

Now  consider  the  stress  intensity  factors  on  the  interface. 

From  Eqn  (3.55),  we  obtain: 


■*  gn(H„r)  +  (H„r)  =  -ij‘  1 JS®  ds 

Tti  J~l  s  —  r 


P1P2P3  2Ej 


P3  2G^ 


(E.2) 


(E.3) 


Irl  >1 


XXIV 


F3°(s)e^ 

where  Vl(s>=  (s+W(t-lf 
Again  from  [33],  we  obtain: 

J.  f1  VMds  =  1  -f?(r)  +  F  (r) ) 

7C1 s-r  (r  +  l)^(r-lr 

where  Fn(r)  is  the  principal  part  of  \|/j  (r)  at  infinity  and  is  bounded.  Therefore: 


F,°(r)] 


Multiply  both  sides  of  the  above  equation  by  -i,  then  we  have: 

^  W  - 

Now  if  we  define  the  following  stress  intensity  factors  at  the  interface: 
1  1 


P3  20^ 


P1P2P3  2E* 


=  Urn  (r+ir(r-l)»[^--i-olwCH1j')  -  i,P^-  °ln(Hi,r)  ] 

Ml  P3  20^  V  P1P2P3  2Ex 


then  we  obtain: 


1  1  -k2(±i)  -  i  f^7 1  y*  k,(±i) = f  Vh^i?  (±d 


(E.4) 


P3  iGjjy  '  V  PiP2P3 

In  order  to  derive  the  expression  for  the  strain  energy  release  rate  at  the  interface 
crack  tip,  we  first  express  \j/3  (r)  in  terms  of  kj  and  k2 ,  then  from  (E.4),  we  have 


P4 


(r+mtr-D^e^.l  1 

Vs  -Jl-C  Pj  2G„  \  P1P1P3 


1  — V  V 
— s^-kj] 


2E. 


r  1  1  , 

[— T7r-k 2  *  1, 


P4 


VRT  Ps  2G*y 


P1P2P3 


1  — V  V 

- S^-K-kJ 


2E. 


f— r 

U  + 


£  1  __L_ 

rj  Vl-r2 


irl  <  1  (E.5) 


Noting  also  that: 


XXV 


1  1 

* 


p4 


P3 *  2EX  V  p!p2p3 


c. 


Lim(r+ir“>(r-l)-|i>[- — —  k,  .  , 

P32G„  VPlpjPj  2E 


i.r5ri~.v?v-kJ 


=  [- — —  ka-i  l-^1  I  -i-j —  Irl >  1 

lP32G„  2  )|P,P2P3  2E,  lJU+lJ 

Assuming  T|  «  1,  where  rj  =  r-1,  then  (E.6)  can  be  approximated  by 


(E.6) 


1  1  .  1_  VxyVyx 

- a„,  -  1 - 5L-^ 


P4 


p3  2G,,  *  2E, 


P1P2P3 


„rl  1  , 

=  [ - k9  - 1 

P>2G„  2  1 

1 


0,  1-V„V„ 

_P4 - SL_a.ki] 


P1P2P3  2Ex 


-  ^  r4cos(mlogOj))  +  i  sin(colog(-j))]  (E.7) 

Similarly,  letting  q  =  1-r  and  assuming  that  q  «  1,  then  (E.5)  can  be  expressed  as: 


\|/°(<; )  s  -  -t===[— — - — k2  -  i  l--?4 . Lv.gvg.k  ] 

P3  1 P1P2P3  2Ex  lJ 

X  ^[cosCcologei))  +  i  sin(colog(-|))] 


(E.8) 


Since  \|/3(r)  =  ^^~[(vrv2)  -  i  J^(uru2)] , 
2  dr  Y  ri2 

then 


YsO; )  =  -  7^[(vrv2)  -  \Ff(uru  2)] 
2  dq  V  ri2 


*  r  1  1  i-  P4  1  v*yvy*  i.  i 

VRr  P>2G„  2  iPlp2p,  2E,  lJ 


XXVI 


:  [cos(cdlog(~ ))  +  i  sin(colog(-|))] 


(E.8*) 


The  oscillatory  behavior  in  (E.8*)  can  be  removed  using  the  techniques  described 
in  [40,42].  Considering  the  fact  that  the  stress  intensity  factors  and  the  displacement 
derivatives  just  away  from  the  interface  crack  tips  should  have  definite  values.  Therefore 
we  can  chose  rj  and  q  in  such  a  way, 

colog(— )  =  0  and  colog(-^)  =  0 
2  2 

In  fact,  for  most  material  combinations,  that  above  two  quantities  are  indeed  very  small 
Then  from  eqns  (E.7)  and  (E.8*),  we  have 


axy  — 


_J__k 

4i^\2 


Av  =  (vrv2)  =  2V2 


1  1  1 


aFCT  P3  2Gv 


and 


Au  =  (uru2)  =  2V2 - 


1  i-Wj, 


^Ftrp,  2B,  ' 

Assume  the  crack  front  was  closed  along  the  interface  by  an  amount  Aa 

then  using  the  method  given  in  [44],  we  obtain  the  increment  of  the  strain  energy: 

<*r+Aa  .  .  .  „ 

AE  =  2Jr  [c^Aufr-Aa)  +  a:tyAv(r-Aa)]dr 
2 


=  2 


t  JLi_VV  k2+ 

»-  ATI  1 


*y  Ty»  1,2  1  ^ 


■JRT  P2  2E 


=  2- 


2  1  l-v»v 


Ps  20^ 
2  1  1 


k2lf 


•x+a«  Vl+ Aa-r 

Vr-1 


dr 


k?+  — — —  k*]-Aa 
1  P3  2G,y  zJ2 


V/RT  P2  2Ex 
Therefore,  the  strain  energy  release  rate  will  be: 

AE  2rc  k2+J — i_k2] 

L  /A  TA  1  _  A  “  -* 


Aa  -yJl  —  ^2  P2  2Ex 


p3  2G,, 


(E.9) 


For  Isotropic  materials  we  have: 


xxvu 


r 

* 

—  --22 — ! — —  [  k;+k|] 

ia  p3  2G„ 

ii)  Case  of  transverse  crack  touching  the  interface 
Define  the  stress  intensity  factor 

k(a,)  =  Lim  2Y(x2+H2)‘>'a2y(x2,0) 

from  eqn  (3.29b),  we  have 

a2y(Xl,0)  =  /i2Yl4-Ey . --  f  K21s(x2,t)(t>1(t)dt  +  a°y(x2,0) 

where  a2y(x2,0)  is  bounded  frmction  and 

am,  F’(t)  _  Fl(t)giY 

9lW  (Hj2  -t2)7  (t-Hj^Ct+Hj)7 


(E.9*) 


(E.10) 


(E.11) 


Again  define: 
then  we  have: 


tp(z)= 


I  p  Ml 

^  J-H,  t-Z 


dt 


-M _ _ - +  cp  (z) 

w  J  (211,  )T  sin  rcy(z+ Hj  )T  (2H1)'sinicy(z-HI)’' 

Substituting  (E.12)  back  into  (E.  1 1),  and  then  from  the  definition  of  (E.10), 


(E.12) 


k(aj)  = 


(H^E^d)  k!/PS 

101  Wilwj/Ps]7 


+  X 


W 


J 


103  [|w3Jwr|/p5]T  101[|w3||w;|/p5] 


wJ/Ps 


-] 


(E.13) 


L1  H  I  -  J  L1 

The  stress  intensity  factors  at  other  crack  tips  and  the  strain  energy  release  rate  at 
the  interface  crack  tip  remain  the  same. 


iii)  H-shaped  cracks 

The  derivation  of  stress  intensity  factors  and  of  the  strain  energy  release  rates  at 
the  interface  crack  tip  remain  the  same  as  those  given  in  section  i). 


